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ABSTRACT 

The  continuing  development  of  the  free  electron  laser  (FEL)  as  a  powerful  and 
versatile  source  of  coherent  radiation  steadily  drives  toward  the  goal  of  high  efficiency 
and  broad  tunability  at  shorter  wavelengths.  New  experiments  provide  significant  data 
and  insight  for  analysis  by  theoreticians  and  experimentalists.  Two  important  areas  of 
study  are  short  electron  pulse  effects,  and  the  dynamics  of  optical  mode  distortion  by 
intense  electron  beam  currents. 

The  initial  part  of  this  thesis  examines  one  aspect  of  the  projected  task  of  FEL 
application  as  a  military  weapon.  The  advantages  of  the  FEL  over  other  directed 
energy  sources  are  detailed,  as  well  as  the  challenge  presented  by  the  effects  of  the 
marine  atmosphere  to  high  energy  laser  propagation. 

The  remainder  of  this  thesis  examines  several  effects  of  long  wavelength  FELs. 
Chapter  IV  examines  the  proposed  parameters  of  the  CEBAF  IR  FEL.  and  the  analysis 
leads  to  predictions  describing  system  performance.  Chapter  V  examines  the  effects 
of  single  pass  optical  mode  distortion  for  FELs  with  narrow  electron  beams.  Single- 
mode  theory  states  that  gain  is  proportional  to  the  product  of  electron  beam  current 
and  filling  factor,  but  three  dimensional  simulations  show  that  gain  is  a  function  of 
electron  beam  filling  factor  alone.  Also  examined  is  a  phenomenon  of  destructive 
interference  of  light  in  the  FEL  undulator. 

Chapter  VI  extends  the  analysis  of  the  relationship  of  gain  and  beam  size  to 
include  multiple  passes  of  laser  light  through  the  laser  resonator.  This  affirms  the 
general  gain  relationship,  where  gain  is  a  function  of  electron  beam  filling  factor,  and 
also  further  explores  the  phenomenon  of  destructive  interference  within  the  optical 
mode. 
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I.  INTRODUCTION 

A  free  electron  laser  (FEL)  produces  light  through  the  stimulated  emission 
from  a  high  energy  electron  beam  in  a  periodic  magnetic  field.  Proposed  by  J. 
M.  J.  Madey  in  1970  [1],  the  FEL  has  great  potential  for  development  as  a 
highly  efficient,  powerful  coherent  light  source  with  medical,  industrial,  and 
military  applications. 

Theoretical  models  and  technological  demonstrations  are  extending  the 
limits  of  FEL  performance,  with  high  power  and  efficiency  possible  within  a  few 
years.  The  FEL  is  a  significant  tool  because  of  its  unique  operating 
characteristics.  Among  these  are  the  use  of  electricity  as  the  power  source, 
order  of  magnitude  tunability  of  output  wavelength,  and  high  "wallplug" 
efficiency.  With  the  likely  development  of  FEL  power  into  the  tens  of  megawatt 
range,  the  potential  exists  for  military  applications.  Ground  based  and  space 
based  options  have  been  studied  by  the  Strategic  Defense  Initiative  Office 
(SDIO).  The  option  of  arming  a  naval  vessel  with  a  high  energy  laser  as  an  air 
defense  weapon  deserves  study,  and  is  discussed  here  in  Chapter  II. 

Chapter  III  provides  a  foundation  of  FEL  theory  to  enable  detailed 
descriptions  of  FEL  performance  in  subsequent  chapters.  Chapter  IV  examines 
the  infrared  (IR)  FEL  proposed  for  the  Continuous  Electron  Beam  Accelerator 
Facility  (CEBAF),  of  Newport  News,  Va.  The  planned  accelerator  and  FEL 
parameters  are  analyzed  to  assist  design  performance  optimization.  The  tuning 
range  of  the  CEBAF  IR  FEL  is  analyzed  to  determine  the  influence  of  short 
electron  pulse  effects  on  FEL  gain  and  power. 


Chapter  V  treats  optical  diffraction,  and  the  phenomenon  of  optical  mode 
distortion  caused  by  intense  electron  beams.  Single-mode  gain  theory,  which 
neglects  the  effects  of  diffraction  and  mode  distortion,  is  shown  to  be 
inadequate.  Using  three-dimensional  simulations  of  the  optical  mode  wave 
front,  a  relationship  between  single-pass  gain  and  electron  beam  size, 
accounting  for  diffraction  and  distortion,  is  found.  A  curve  is  generated  which  is 
generally  applicable  to  all  FELs  with  small  electron  beams.  Chapter  VI  explores 
the  new  gain  relationship  and  predicted  destnjctive  interference  phenomena  in 
a  multiple-pass  simulation.  A  similar  gain  curve  is  generated. 


II.  SHIPBOARD  FREE  ELECTRON  LASER 

1.  INTRODUCTION 

The  current  state-of-the-art  anti-ship  missile  (ASM)  threat  posed  to  the  U.S. 
Navy  consists  of  high  speed,  maneuverable,  robust  and  electronically 
sophisticated  weapons.  Current  Fleet  defensive  weapons  and  tactics  center  on 
the  highly  complex  Aegis  Combat  Systems  Suite,  which  integrates  sensors,  fire 
control,  and  decision  aids  with  surface-to-air  missiles  (SAMs).  The  Achilles' 
Heel  of  this  system  is  the  limited  number  of  SAMs  and  their  performance 
limitations.  High  Energy  Laser  (HEL)  design  concepts  have  been  studied  for 
incorporation  into  combat  ships  for  air  defense.  The  Free  Electron  Laser  (FEL), 
due  to  its  design  and  performance  characteristics,  would  be  the  optimal  choice 
for  a  future  ship  based  HEL  weapon. 

2.THREAT  ENVIRONMENT 

Combating  the  high  speed  and  possible  high  density  of  anti-ship  guided 
missiles  is  among  the  highest  priority  missions  of  tactical  naval  planning  for  the 
next  generation  of  shipboard  weapon  systems.  ASMs,  capable  of  launch  from 
submerged  submarines,  tactical  aircraft,  and  small  patrol  boats  are  cheap, 
reliable  and  difficult  to  defeat  with  all  but  the  most  complex  integrated  systems. 
In  worst-case  wargaming  scenarios,  which  include  numerous  ASMs 
simultaneously  converging  from  many  directions,  even  the  most  advanced  air 
defense  system,  such  as  an  Aegis-equipped  guided  missile  cmiser,  can  be 
overwhelmed  through  exhaustion  of  the  SAMs  carried.    This  is  known  as  the 


"magazine  depth." 

Time  is  also  considered  a  highly  limiting  factor  for  anti-ship  missile 
defense.  With  present  systems,  once  a  target  is  identified,  evaluated  as  a 
threat  and  the  order  is  given  to  engage,  there  is  a  time  lag  due  to  SAM  flight  to 
intercept  the  target  lasting  from  ten  seconds  for  a  close  range  engagement  to 
two  minutes  for  a  longer  range  target.  If  the  target  is  not  destroyed,  this 
process  must  be  repeated,  further  wasting  time  and  valuable  defensive  missiles. 
The  Aegis  system  can  successfully  track  and  engage  several  incoming  missiles 
simultaneously,  but  is  ultimately  limited  by  its  magazine  depth. 

3.  HEL  OPTIONS 

HELs,  having  been  the  subject  of  extensive  research  since  the  initiation  of 
the  Strategic  Defense  Initiative  (SDI)  program  in  1983,  are  usually  categorized 
by  the  fuel  sources  that  drive  them.  These  include  X-Ray  Lasers,  Excimer  or 
Chemical  Lasers,  and  Free  Electron  Lasers. 

X-Ray  Lasers,  despite  characteristics  desirable  for  exoatmospheric  SDI 
missions,  are  not  a  feasible  option  for  a  shipboard  HEL.  The  wavelength  of 
radiation  produced  by  these  lasers,  in  the  X-Ray  region,  does  not  propagate 
effectively  in  the  atmosphere,  rendering  it  somewhat  useless  as  a  shipboard 
weapon.  More  importantly,  the  laser  is  powered  by  the  detonation  of  a  nuclear 
weapon,  which  understandably  obviates  further  discussion  of  its  adoption  by  the 
fleet. 

Chemical  lasers,  having  been  tested  successfully  as  HELs,  are  a  much 
more  attractive  option  for  a  shipboard  weapon.  The  MIRACL  laser  system, 
employing  a  deuterium  fluoride  (DF)  laser,  successfully  destroyed  a  target 


drone  aircraft  at  White  Sands  Missile  Range.  This  laser  emits  light  in  a 
favorable  wavelength  for  atmospheric  use  (  3.7-4.0^, )  [2],  and  has  been  studied 
for  possible  shipboard  use.  The  main  obstacle  for  use  of  chemical  lasers  is  the 
storage  and  combustive  reaction  of  the  highly  volatile  agents  needed  to  drive 
the  laser.  The  exhaust  gases  of  chemical  lasers  are  highly  toxic  and  would 
require  a  complex  air  filtration  and  cleanup  system  to  enable  habitability  for  the 
crew  of  a  vessel  equipped  with  a  chemical  laser. 

Chemical  lasers  are  driven  through  the  exothermic  reactions  of  highly 
volatile  oxidizer  and  fuel  streams  which  are  mixed  inside  the  optical  cavity.  The 
reactive  nature  which  enables  these  agents  to  be  useful  is  itself  a  strong  reason 
for  the  rejection  of  chemical  lasers  as  shipboard  weapons.  The  transport  and 
storage  of  highly  pressurized,  toxic,  and  reactive  chemicals  to  ships  at  sea 
would  require  designs  of  replenishment  ships  (tankers  and  supply  ships)  that 
today  do  not  exist.  The  potential  for  disaster  once  the  chemicals  were  aboard 
the  laser  ship  would  necessitate  detailed  and  expensive  safety  programs, 
further  complicating  the  process.  Ultimately,  the  chemically  driven  laser  system 
would  fall  victim  to  the  same  limitations  of  the  missile  systems  it  would  seek  to 
replace.  The  chemical  laser  is  driven  by  an  irreversible  chemical  reaction  that 
would  cease  once  the  fuel  and  oxidizer  were  exhausted.  This  is  a  solid 
magazine  depth  that  would  severely  limit  the  duration  of  engagement  and 
require  complex  reloading  procedures.  This  is  a  major  flaw  for  a  plan  to  field  a 
chemical  HEL  aboard  a  ship.  A  principle  advantage  of  a  ship  as  a  weapons 
platform  includes  long  endurance  on  station,  and  it  would  be  wasteful  of  this 
endurance  potential  to  place  a  short-legged  weapon  aboard  as  a  principle 
means  of  air  defense  for  a  battle  group  or  geographic  region. 


A  simpler  and  more  logical  system  for  deployment  of  HELs  for  near  and 
long  term  use  would  involve  development  of  an  aircraft-based  chemical  HEL, 
followed  later  by  the  development  of  a  shipboard  FEL  (SFEL)  when  the 
necessary  technology  matures.  Aircraft  lack  the  endurance  of  ships,  needing  to 
land  frequently  for  fuel,  maintenance  and  crew  rest,  so  that  the  engage-and- 
reload  requirements  for  a  chemical  HEL  would  not  represent  as  severe  a 
mission  limiting  factor  as  would  be  the  case  for  a  shipboard  chemical  HEL.  The 
lessons  learned  through  development  of  the  aircraft  based  HEL  would  later 
benefit  the  deployment  of  the  SFEL 

The  Free  Electron  Laser  has  several  advantages  over  other  potential  HEL 
designs.  The  principle  advantage  of  this  system  is  the  wide  range  of 
wavelengths  over  which  the  FEL  can  be  tuned.  This  allows  the  FEL  to  be 
tuned  to  exploit  highly  variable  atmospheric  conditions  and  to  be  more  highly 
resistant  to  countermeasures.  An  additional  advantage,  crucial  when 
considering  distribution  of  large  quantities  of  electricity  aboard  ships,  is  that 
scaling  up  to  allow  the  addition  of  an  FEL  would  be  simply  a  matter  of  adding 
more  generators  and  power  conditioning  equipment.  Recent  technological 
advances  in  capacitor  storage  banks  would  allow  large  quantities  of  electrical 
energy  to  be  held  in  reserve  for  immediate  use,  providing  a  quick  reaction 
capability  for  the  SFEL.  This  is  an  easily  manageable  task,  and  would  fall  in 
well  with  the  navy's  future  plans  for  an  integrated  electric  drive  combatant 
vessel.  An  especially  attractive  factor  of  FEL  design  is  the  generally  accepted 
"wallplug"  efficiency  goal  of  40%.  This  is  would  be  made  possible  through  use 
of  an  energy  recovery  system,  which  would  also  enable  reduction  of  the 
necessary  mass  of  radiation  shielding.  The  limiting  factor  for  the  operation  of 
the  FEL  would  then  be  the  common  fuel  supply  used  for  propulsion  and 


electricity  of  the  ship,  which  is  already  supported  by  a  well  established  system. 

4.  ATMOSPHERIC  PROPAGATION  EFFECTS 

The  marine  atmospheric  environment  within  which  an  SFEL  would  operate 
is  a  highly  variable,  dynamic  system  with  numerous  linear  and  nonlinear  effects 
that  have  profound  influence  on  the  propagation  of  high-power  laser  beams.  In 
addition  to  the  effect  of  optical  diffraction  within  the  atmosphere,  there  exist 
atmospheric  attenuation  effects.  The  mechanisms  of  laser  beam  photon  loss 
through  atmospheric  attenuation  by  the  constituents  suspended  in  the 
atmosphere  can  be  broken  down  into  several  distinct  effects:  the  linear 
absorption  and  scattering  of  the  molecular  and  aerosol  elements  in  the 
atmosphere,  the  random  wander,  spreading,  and  distortion  of  the  laser  beam 
due  to  the  effects  of  turbulence  within  the  atmosphere,  and  the  non-linear 
defocusing  effects  of  thermal  blooming  caused  by  the  absorption  of  small 
amounts  of  beam  power  by  the  atmosphere.  These  effects  are  dependent  on 
atmospheric  temperature,  clarity,  humidity,  wind  velocity,  and  beam 
characteristics  such  as  wavelength,  temporal  operating  mode  and  beam 
diameter. 

Accounting  solely  for  linear  atmospheric  effects,  the  mean  peak  irradiance 
of  a  propagating  laser  beam  is  Ip  =  [P/{nw^)]  ex^i-ajz),  where  P  is  the 
transmitted  laser  power,  z  is  the  range,  and  ar  is  the  atmospheric  extinction  or 
attenuation  coefficient.  The  term  w  is  the  root  mean  square  (rms)  beam  radius 
and  is  a  function  of  many  factors,  w  =  yjw^  +  wf  +  wf,  where  w^  is  a  factor  of 
diffraction  and  beam  quality,  wy  accounts  for  beam  jitter,  and  w,  is  a  correction 
for  atmospheric  turtDulence  [3]. 


The       diffraction       and       beam       quality       factor       correction       is 


v/a  =  ^'^^z^ik^vj'^  +  w\\  -  ziRf.    Here  p  is  a  ratio  of  actual  focal  spot  radius  to 
the  diffraction  limited  spot  radius,  or 

R  _    1/e  Radius  of  actual  focal  spot  _     ^d  _  ^/-x,^  o  -• 

^,        P  "  1/e  Radius  of  Dlff.  Limited  spot      w,^^  '       '^^  "     ^  ' 

Range  R  can  be  varied  from  R  =  z  for  a  beam  focussed  to  the  point  at 
distance  z,  to  /?  =  ~  for  a  collimated  beam.  Beam  jitter  effects  are  corrected  for 
by  the  addition  of  the  mean  square  radial  displacement  of  the  focal  spot, 
Wj  =  V2<0j>?,  at  range  z  for  a  varying  single  axis  jitter  angle,  <0^.  This  angle 
is  isotropic  for  x  and  y  axes,  so  that  <e|>  =  <e^>. 

The  attenuation  and  extinction  term,  a^  =  aa^s  +  otscaf.  represents  the 
effects  of  absorption  and  scattering  of  molecular  and  aerosol  components  of  the 
atmosphere.  Scattering  is  a  directional  redistribution  of  the  beam  incident  on 
the  various  molecules  composing  the  atmosphere.  This  does  not  reduce  the 
energy  of  the  beam  but  does  reduce  the  intensity.  Absorption  is  a  very 
complicated  and  irregular  function  of  beam  wavelength,  and  is  a  critical  factor  in 
the  effects  of  thermal  blooming.  Scattering  and  absorption  by  aerosol 
components  is  generally  less  dependent  on  wavelength  than  molecular 
absorption  and  scattering.  The  slight  wavelength  dependence  of  the  absorption 
and  scattering  coefficients  is  a  result  of  the  combined  effects  of  the  aerosol 
complex  refractive  index  and  the  size  distribution  of  the  aerosol  particles. 
Aerosol  effects  become  somewhat  significant  at  shorter  wavelengths 
(><3|i  m)[3].  Fig  2-1  shows  the  highly  variable  nature  of  the  aerosol  and 
molecular  effects. 
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Figure    2-1:    Calculated    aerosol    absorption    and    extinction    by   the 
atmosphere  in  23km  visibility  vs.  wavelength  [3]. 


Figure  2-1  is  generated  through  a  model  of  the  continental  atmosphere. 
The  marine  environment  would  include  higher  levels  of  humidity  and  its  effects. 
Included  on  the  graph  are  range  bars  indicating  the  wavelengths  of  three 
common  lasers  and  the  variation  of  the  attenuation  coefficients  for  each.  Of 
note  is  the  dip  in  aerosol  absorption  effects  for  X  =  3.8p.m,  which  closely 
corresponds  to  the  wavelength  of  the  DP  laser.  This  wavelength  is  near  an 
optimum  wavelength  selected  for  the  study  of  FEL  placement  on  a  naval  vessel 
[4]. 

Atmospheric  transmittance,  accounting  for  molecular  and  aerosol 
absorption  and  scattering  effects,  is  a  ratio  of  initial  irradiance,  /q.  to  that  at  a 
range  z  : 


r  =  igi  =  exp{-J>z}. 


2-2 


Atmospheric  turbulence  arises  from  convective  air  currents  and  eddies 
caused  by  air-surface  temperature  gradients,  which  result  in  slight  fluctuations 
in  the  index  of  refraction.  For  beam  propagation,  the  eddies  form  pockets  of 
disturbed  air,  known  as  "turbules",  which  vary  in  size  relative  to  the  beam 
diameter.  This  non-uniformity  results  in  the  various  effects  of  beam  wandering 
off  the  target,  intensity  fluctuations  across  the  beam  profile,  and  spreading  of 
the  beam  spot  size.  Turbulent  beam  spread  and  wander  effects  are  accounted 
for  by  the  term  w,  =  2z/kLo  where  z  is  distance,  k  is  the  wave  number,  and 
Lo  =  {0.545k^ C^z)~^'^  is  the  lateral  coherence  length  of  a  spherical  wave.  C/J  is 
the  refractive  index  structure  constant.  Figure  2-2  shows  wavelength 
dependence  of  the  effects  of  beam  spreading  and  wander  on  average 
irradiance,  l{X),  normalized  by  the  irradiance  at  X=  tOp-m  [3].  At  long  ranges, 
turbulent  beam  spreading  becomes  wavelength  independent,  but  this  is  an 
exceptional  regime.  In  the  usual  case,  the  turbulence  coefficient 
w,  =  2.0^(k-^'^C^^z^^).  Though  the  dependence  of  turbulent  beam  spread  on 
X~^'^  is  relatively  weak,  this  relationship  implies  that  there  is  an  optimum 
wavelength  for  propagation  for  given  values  of  turbulence  level,  range,  and 
aperture  diameter  (diffraction  being  proportional  to  X).  Neglecting  jitter  effects, 
the  mean  peak  intensity  is  given  by 

/  oc  J ^ .  2-3 

w^     1  -I-  Wt^/W^ 

Assuming  a  reference  wavelength  Aq,  the  ratio  of  the  irradiances 

2-4 

J     1  -h>A(Ao/X)"^=' 

where  A  =  16.l27c^vv^(C/vzA^^)"^'^/p^.    For  wavelengths  shorter  than  the  critical 


/(>>)    =l^\  1  +A 

/(Ao)       [Xj     1+>\(Ao/X)^2/5 
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wavelength,  defined  by  X^  =  {Aisf^'^h^,  the  effects  of  turbulent  spreading 
dominate  the  diffraction  limited  focusing,  and  beam  intensity  decreases  as  X^^ 
[3]. 
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Figure  2-2:  Dependence  of  turbulence  effects  on  wavelength  [3]. 


Thermal  blooming  or  defocusing  is  a  time  dependent  effect  of  absorption  of 
a  small  amount  of  beam  energy  by  the  atmosphere  along  the  beam  path.  The 
resultant  rapid  heating  lowers  the  air  density  and  also  the  refractive  index, 
causing  a  nonlinear  negative  lens  effect  of  divergence  on  the  beam.  Varying 
with  atmospheric  conditions,  defocusing  can  occur  in  as  little  as  0.2s.    This 
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length  of  time  is  crucial  because  it  is  close  to  a  typical  value  of  the  "dwell  time" 
necessary  for  a  laser  beam  to  remain  focussed  in  order  to  damage  or  destroy  a 
target.  The  spreading  effect  can  dominate  diffraction  limited  focusing.  This  can 
limit  the  irradiance  on  a  target  independent  of  the  available  power  of  the  beam. 
Accounting  for  variable  atmospheric  velocity  due  to  wind  and  the  slew  rate  of  a 
beam  tracking  a  target,  there  will  be  a  temperature  gradient  across  the  face  of 
the  beam,  with  the  "cooler"  air  on  the  "upwind"  side  of  the  beam  face,  in  the 
direction  of  the  cumulative  wind/slew  velocity.  The  temperature  gradient  will 
result  in  a  varying  refractive  index  across  the  beam  face,  leading  to  an 
asymmetric  distortion  and  growth  of  the  beam  cross  section.  With  the  "upwind" 
side  of  the  beam  face  remaining  at  a  lower  average  temperature,  this  side  will 
tend  to  have  a  higher  air  density  and  thus  a  higher  index  of  refraction.  This 
causes  the  path  of  the  defocused  beam  to  bend  "upwind"  asymmetrically. 

Peak  irradiance  on  a  target  accounting  for  thermal  bloom  effects  must  be 
scaled  by  a  factor  of  thermal  distortion,  N,  which  is  a  function  of  focusing  of 
beam  radii,  the  effects  of  a  non-uniform  velocity  profile,  and  atmospheric 
attenuation,  aj.  The  specific  heat,  Cp,  of  the  air,  the  range,  z,  and  the  power, 
P,  of  the  laser  are  also  accounted  for.  Essentially,  N  «  ajPoZivD^,  where  v  is 
a  wind  velocity  term,  a©  is  the  collimated  beam  Me  radius,  and  D  =  2V2ao  is  the 
aperture  diameter.  An  empirical  relationship  that  works  well  in  modeling 
systems  evolves  as  /re/ =  Wom/^peaA  =  1/(1  +  0.06A/^).  Since  peak  irradiance  is 
proportional  to  the  product  of  power  and  the  empirical  correction,  Ipgg,^  «  Px/^g/, 
the  peak  irradiance  on  target  will  increase  until  /^e/  decreases  at  a  rate 
proportional  to  1/P  or  greater.  This  indicates  that  for  a  given  beam  size,  focal 
range  and  atmospheric  conditions,  there  exists  a  peak  power,  P^.,  at  which 
irradiance  on  target  is  maximum. 
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Using  this  empirical  relationship,  Fig  2-3  shows  the  effects  of  thermal 
blooming  in  limiting  peak  irradiance  and  maximum  power  for  a  propagating 
beam. 


P/Pc 


Figure  2-3:  The  normalized  peak  irradiance  vs.  the  normalized  power  PiPq, 
based  on  the  empirical  relationship  of  Ibjoom^lpeak  [3]. 


Unlike  other  atmospheric  effects,  increased  power  in  the  laser  beam  will 
not  result  in  higher  fluence  on  the  target,  but  will  only  serve  to  further  aggravate 
the  blooming  effect.  This  is  a  nonlinearity  unlike  any  other  atmospheric  effect. 
All  other  phenomena  result  in  increased  target  irradiation  with  higher  power. 
Peak  in-adiance  accounting  for  linear  and  nonlinear  effects  at  range  R  =  z  \s 

P       exp(-a72) 


/p  = 


na^    1  +  0.0625/V2  ' 


2-5 
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The    linear   effects   of   diffraction,    beam   quality,   jitter,   and   turbulence   are 
accounted  for  by  the  rms  beam  radius 

w2  =  ^^yf  +  2<e2>z2  +  AiC^zf^X-^^z^  2-6 

where  D  -  2V2ao  is  the  aperture  diameter,  and  the  linear  attenuation  effects  are 
covered  by  07  =  03^,5  +  a^cax- 

The  effects  of  thermal  blooming  are  time  dependent,  and  so  there  is  a 
difference  between  the  propagation  effects  for  continuous  wave  (cw)  and  pulsed 
beams.  For  a  pulsed  beam  of  diameter  D,  if  the  pulse  repetition  frequency 
(PRF)  is  less  than  the  time  needed  for  the  cross  wind,  »^o.  to  clear  the  tube  of 
atmosphere  heated  by  the  previous  pulse,  or  PRF<i^o/^.  there  will  be  no 
thermal  bloom  effects  from  overlapping  pulse  beam  effects.  For  time-scales 
outside  of  this  regime,  however,  both  pulsed  and  cw  beams  will  experience 
bloom  effects  due  to  continuous  heating.  Through  complex  studies  of  various 
PRFs  and  beam  face  irradiance  levels  for  pulsed  beams,  it  has  been  found  that 
pulsed  beams  are  less  affected  by  atmospherics  than  cw  beams  at  long  ranges 
[3]. 

Marine  atmospheric  propagation  effects  are  difficult  to  predict  on  any  useful 
time  scale.  The  few  controllable  parameters  include  light  wavelength  and 
power  level  of  the  laser.  Chemical  HEL  designs  studied  for  naval  application 
have  all  had  basic  flaws  of  limited  fuel  capacity  and  fixed  wavelength.  Free 
electron  lasers,  though  less  mature  in  design,  have  the  valuable  characteristics 
of  tunable  wavelength  and  manageable  fuel  and  power  requirements.  The 
minimizing  of  atmospheric  effects  on  laser  beam  propagation  will  need  to  be 
carried  out  not  only  by  the  selection  of  laser  parameters,  but  also  in  creative 
beam  optics.    A  system  capable  of  atmospheric  propagation  diagnostics  and 
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adaptation  of  deformable  optics  would  greatly  enhance  the  performance  of  any 
high-energy-laser  weapon.  Atmospheric  diagnostics,  hitherto  the  realm  of 
astronomical  studies  for  adaptive  telescope  optics,  should  be  developed  for 
inclusion  into  the  control  system  of  any  SHEL  design.  In  astronomy,  optical 
system  diagnostics  are  performed  by  observation  of  a  "guide  star"  of  known 
characteristics.  For  tactical  military  purposes,  a  low  power  laser  carried  by 
drone  aircraft  or  other  suitable  platform  could  be  aimed  to  illuminate  diagnostic 
optics  to  allow  analysis  of  the  prevailing  atmospheric  conditions  and  subsequent 
optimization  of  the  beam  parameters  and  adaptive  optics'  adjustment  [2,5]. 

A  more  exotic  solution  may  be  found  by  stationing  a  relay  mirror  system 
high  above  the  FEL  vessel  aboard  a  long-endurance  fixed-wing  remotely-piloted 
vehicle  (RPV)  or  balloon  aircraft.  A  system  centered  on  this  concept,  named 
Thunderball,  was  presented  to  the  Physics  department  of  the  Naval 
Postgraduate  School  by  Lieutenant  Colonel  Ed  Pogue,  USA,  former  Deputy 
Director  of  Directed  Energy  for  the  Strategic  Defense  Initiative  Office.  This 
concept  had  been  studied  by  the  Advanced  Technology  Group  of  W.  J.  Schafer 
Associates  [6].  The  chief  advantages  of  this  system,  whose  aircraft  would  be 
stationed  at  altitudes  approaching  eighty  thousand  feet,  would  be  the 
minimization  of  distance  traveled  through  the  dense  lower  atmosphere.  This 
would  help  minimize  the  most  harmful  effects  of  the  atmosphere.  Additionally  a 
beam  of  greater  diameter,  and  thus  lower  power  density  across  the  beam  face, 
could  be  used.  This  would  lessen  the  damaging  effects  of  atmospheric 
influences.  Once  the  beam  reached  the  optical  system  aboard  the  aircraft,  a 
beam  cleanup  and  refocusing  system  could  narrow  and  redirect  the  beam 
toward  a  target  far  beyond  the  horizon  of  the  FEL  vessel  below.  This  system 
would  be  well  suited  for  a  defensive  mission  against  incoming  ballistic  missiles. 
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III.   FREE  ELECTRON  LASER  THEORY 

A.  BASIC  FREE  ELECTRON  LASER  PHYSICS 

Free  electron  laser  systems  consist  of  four  basic  components,  the  electron 
beam  accelerator,  the  optics  which  include  the  beam  transport  system  and  the 
PEL  undulator  and  resonator.  The  accelerator  produces  a  stream  of  relativistic 
electrons  that  is  injected  into  the  periodic  magnetic  field  of  the  FEL  undulator, 
which  "wiggles"  the  electrons  along  sinusoidal  trajectories.  The  periodic 
transverse  acceleration  of  the  highly  relativistic  electrons  results  in  classical 
synchrotron  radiation  into  a  forward  cone  along  the  beam  path.  The  resonator 
cavity  consists  of  two  curved  mirrors,  one  is  fully  reflecting  and  one  is  partially 
reflecting,  placed  beyond  the  ends  of  the  axis  of  the  undulator.  The  mirrors  are 
designed  to  sustain  a  transverse  Gaussian  mode  shape.  The  total  resonator 
length  is  typically  S  =  5  -  20m.  The  FEL  oscillator  starts  from  spontaneous 
emission,  some  of  which  is  saved  in  the  resonating  cavity  resulting  in  an 
increase  in  optical  power.  Repeated  coupling  of  the  electrons  with  the  optical 
and  undulator  fields  results  in  the  bunching  of  the  electrons.  This  amplifies  the 
amplitude  of  the  wavefront,  resulting  in  coherent  radiation  and  lasing.  Undulator 
gain  is  the  fractional  change  in  optical  power  per  pass,  G  =  AP/P.  Typical 
values  for  G  =  0.10  for  a  low  gain  FEL  and  G  =  10^  for  a  high  gain  FEL. 

The  bunching  effect  is  cmcial  to  FEL  operation,  and  is  highly  dependent  on 
the  quality  of  the  electron  beam  produced  by  the  electron  accelerator. 
Specifically,  the  energy  spread  of  the  electrons  as  well  as  the  angular  spread  of 
the  beam  must  be  reasonably  low  for  FEL  success,  the  energy  spread  among 
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the  electrons  in  the  beam  results  in  differences  of  velocity  within  the  undulator, 
leading  to  an  axial  smearing  of  the  electron  bunch  entering  the  undulator.  If  the 
velocity  differences  spreads  the  electrons  beyond  the  optical  wavelength,  X,  the 
FEL  gain  would  be  seriously  degraded.  Variations  of  the  electron  injection 
angle,  or  angular  spread,  results  in  non-axial  components  of  bunch  velocity. 
Reducing  the  axial  velocity  will  also  serve  to  degrade  gain  in  much  the  same 
way  as  an  energy  spread. 

Emittance  and  beam  spreading,  collectively  known  as  "beam  quality," 
influence  bunching  effects.  Emittance,  e  =  re,  is  the  product  of  the  root-mean- 
square  (rms)  values  of  radial  and  angular  spreads  of  the  electron  beam.  In  the 
tuning  of  a  given  electron  transport  system,  both  of  these  parameters  may  be 
changed,  but  their  product,  e  =  r  8,  remains  constant.  The  normalized 
emittance,  £„  =  ye,  is  the  product  of  the  electron  beam  Lorentz  factor,  y,  and  the 
emittance  e.  The  normalized  emittance  tends  to  remain  constant  even  when 
the  electron  beam  is  accelerated  to  higher  energies.  Normalized  omittances 
typically  vary  from  a  few  tens  to  a  hundred  mm-mrad. 

The  accelerator  produces  relativistic  bunches  of  electrons  with  peak 
currents  from  1  A  to  IkA,  and  bunch  lengths  of  Ig  =  lO'^cm.  The  relativistic 
electron  beam  will  typically  have  a  Lorentz  factor,  y,  from  ten  to  the  low 
hundreds.  Typical  electron  beam  radii  are  r^  ~  0.10cm.  These  result  in  electron 
densities  of  p  =  lO^^cm"^  to  lO^^cm"^.  With  Lorentz  factors  y  varying  from  ten  a 
few  hundred,  the  beam  energies,  (Y-1)mc^,  range  from  several  tens  of  to 
hundreds  of  MeV,  where  m  is  the  electron  mass,  c  is  the  speed  of  light,  and 
mc^  =  0.51 1  MeV  is  the  electron  rest  mass. 
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The  undulator  consists  of  a  series  of  alternating  magnetic  pole  pieces 
forming  an  intense  magnetic  field  that  accelerates  the  relativistic  electron  beam 
into  periodic  transverse  motion.  Within  the  undulator,  electron  bunches  undergo 
one  cycle  along  an  axial  length  known  as  an  undulator  wavelength,  Xo,  which  is 
typically  five  to  ten  cm.  The  number  of  undulator  wavelengths,  N,  varies  from 
twenty  five  to  250,  which  results  in  undulator  length  L  =  A/Xo  from  one  to  ten 
meters.  The  undulator  magnetic  field  strength,  B,  ranges  from  3kG  to  5kG. 
Mean  undulator  magnetic  field  strength,  E,  is  rms  field  averaged  over  the 
undulator  period  Xq.  For  helically  (circularly)  polarized  undulators,  B  =  B,  due  to 
the  constant  acceleration  of  the  electrons  in  a  helical  path  along  the  axis  of  the 
undulator.  Helical  undulators  produce  circularly-polarized  radiation.  Magnetic 
fields  for  linearly-polarized  undulators  result  in  non-constant  acceleration  due  to 
the  sinusoidal  electron  path,  so  that  B  =  S/V2.  Linearly-polarized  undulators 
produce  linearly-polarized  radiation. 

These  parameters,  with  the  electron  charge  magnitude,  e,  combine  as  the 
"undulator  parameter"  [7], 

eBXo 


K  = 


2Kmc^ 


The  undulator  parameter,  K  varies  from  0.1  to  3.  K  may  be  varied  by  altering 
the  separation  distance  between  pole  faces  along  the  undulator  in  order  to 
change  B.  The  optical  radiation  wavelength,  X,  produced  by  the  FEL  is  given 
by  [7] 

The  optical  wavelength  is  a  function  of  the  undulator  magnetic  field  and  the 
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electron  energy  from  the  accelerator,  as  represented  by  the  Lorentz  factor  y. 
This  relation  reflects  the  tunability  of  the  FEL  over  a  very  broad  range  of 
wavelengths.  Experiments  have  a  x10  range  of  tunability,  and  the  FEL  can 
operate  from  the  mm  to  the  X-Ray  wavelength  range. 

Short  optical  wavelengths  are  a  result  of  the  Doppler  shift  effect  on  the 
undulator  and  optical  radiation  fields  as  "seen"  by  the  highly  relativistic 
electrons.  The  electrons  "see"  a  Lorentz  contracted  undulator  wavelength, 
V  =  ^  /y-  The  electrons  emit  light  at  a  longer  Doppler  shifted  wavelength 
X  ~  2iK.  With  the  similarity  of  the  Doppler  shifted  light  and  the  Lorentz 
contracted  undulator  wavelength,  X  ~  V.  the  relationship  of  the  radiation 
wavelength  to  the  undulator  wavelength  becomes  X  =  Xo/2y^-  This  is  the 
mechanism  which  enables  short  wavelengths  to  be  produced. 

The  resonance  condition  is  met  when  one  wavelength  of  light  passes  over 
an  electron  as  the  electron  traverses  through  the  length  of  one  undulator 
wavelength.  As  the  electron  bunch  is  subjected  to  the  resonant  forces  of  the 
combined  undulator  and  optical  fields,  the  electrons  undergo  either  stimulated 
emission  or  absorption,  depending  on  the  phase  relationship  of  the  bunch 
relative  to  the  optical  field.  This  new  light  from  the  electrons  will  interfere  with 
the  light  of  the  optical  mode.  If  the  net  interaction  of  the  radiation  is 
constmctive  interference,  there  exists  a  regime  of  positive  gain,  G  =  AP/P  >  0, 
for  the  undulator. 

Figure  3-1  illustrates  the  interaction  of  an  electron  with  the  optical  mode 
and  magnetic  fields  within  the  undulator.  The  top  diagram  shows  an  electron 
bunch  and  an  optical  pulse  in  the  undulator.  The  middle  diagram  illustrates  the 
details  along  one  undulator  period. 
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Figure  3-1 :  Fundamental  design  and  electron-optical  mode  interactions 
within  the  undulator  of  an  FEL 
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The  electrons  follow  a  sinusoidal  trajectory  along  the  periodic  magnetic 
undulator  field  while  one  wavelength  of  light  passes  over  each  electron.  The 
bottom  diagram  details  the  electron-optical  mode  interaction.  Depending  on  the 
relative  phase  between  the  optical  mode  and  the  electron,  the  electron  will  be 
subject  to  either  a  retarding  or  accelerating  force  that  will  induce  it  to  either  emit 
or  absorb  energy.  Conservation  of  energy  determines  that  the  energy 
transferred  from  the  electrons  must  be  absorbed  by  the  optical  mode  wave. 
This  results  in  optical  amplitude  growth  and  FEL  gain.  Should  the  phase 
relationship  cause  the  electron  to  be  accelerated,  the  force  to  do  so  will  rob  the 
optical  mode  of  energy  and  result  in  decreased  amplitude,  or  negative  FEL 
gain.  In  non-resonant  systems,  gain  can  be  positive  or  negative  in  various 
positions  along  the  undulator  resulting  in  complex  optical  amplitude  profiles. 

B.  ELECTRON  DYNAMICS 

The  highly  relativlstic  electron  beam  is  subject  to  the  combined  forces  of 
the  radiation  of  the  optical  mode  and  the  magnetic  field  of  the  undulator.  The 
equations  of  the  Lorentz  force  exerted  on  the  electrons  are 

where  ]?=  x^/c  is  the  ratio  of  electron  velocity  V  to  the  speed  of  light  c,  ^  is  the 
optical  electric  field,  and  the  net  magnetic  field,  ^  =  5-  +  ^u  is  the  sum  of  the 
magnetic  fields  of  the  optical  mode  and  the  helical  undulator  [7].  Using  the 
combined  undulator  and  optical  fields  in  3-1,  and  integrating  over  time,  the 
transverse  electron  motion  within  the  undulator  may  be  solved.  There  are  a 
total  of  five  equations  in  3-1,  with  only  four  unknown  terms  [x'(r),  y(0].  so  that 
one  of  the  five  equations  may  be  disregarded.   The  relativlstic  electrons  have 
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an  initial  position  Zq,  an  axial  velocity  V^  =^c,  where  p^  =  1  -  (1  +  K^)/-f  =  1, 

and  the  position  along  the  undulator  axis  is  given  by  z{t)  =  ZQ  +  ^^ct  + The 

transverse  velocity  pj^  =  V1  -  Pz  <:  1,  so  the  magnitude  of  the  transverse  motion 
of  the  electrons  is  much  less  than  the  axial  motion.  The  transverse  undulator 
fields  are  much  greater  than  the  transverse  optical  fields,  so  that  the  transverse 
force  equation  of  the  electrons  becomes 

K  =  4(Y  Pi)  =  ^Pz  e[-sin{;co  2).  cos(/co  z).  0].  3-2 

or  mc 

Where  /cq  =  27i/Xo  is  the  undulator  wave  number.  Integrating  this  equation  over 
time,  and  assuming  perfect  beam  injection,  yields  the  transverse  velocity: 

?i  =  —  [  cos(/(o  z),  s\n{ko  z),  0].  3-3 

One  more  integration  will  result  in  the  transverse  position  relationship  for  the 
electrons: 

Xi=YY^  ^[  sin(/co  z).  -<;os(/co  z),  0]. 

It  can  easily  be  seen  here  that  because  K  =  1  and  the  Lorentz  factor  y=  100 
that  the  magnitude  of  the  transverse  motion  is  relatively  small. 

C.  ENERGY  EVOLUTION 

The  second  term  in  the  Lorentz  force  equations  3-1  determines  the 
evolution  of  the  energy.  Using  the  expression  above  for  transverse  velocity  pj^ 
3-3,  the  approximate  electron  position  z  =  Zq  +  ^  ct  +  •  •  • ,  for  Lorentz  factor 
Y»  1,  the  energy  for  a  weak  field  regime  may  be  approximately  calculated.  The 
optical  electric  field  is 
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^r  =  ^  (cos(/cz  -  (Of).  -sin(/cz  -  cof).  0), 

where  the  wave  number  is  given  by  ^  =  27c/X,  and  the  principle  carrier  frequency 
is  given  by  o)  = /cc  =  27cc/X[].  Using  y  =  -e{^  •  ^)/mc,  the  initial  position 
zq  =  zo  +  Po  ct,  and  substituting  the  expression  for  ^;.,  the  energy  term  evolves 
as 

Y  =  -^^cos((  k  +  ko)z  -  (Of)  =  -^(X)S((  k  +  ko)  Zo  +  [{k  +  ko)  %-k]  ct),  3-4 
ymc  ymc 

where  Po  is  the  ratio  of  initial  electron  velocity  vq  and  the  speed  of  light  c  [8]. 
Define  the  initial  phase  as  Co  =  (^o  +  ^)^o  =  27c  Zq/X,  which  relates  the  initial 
electon  z  position  relative  to  the  optical  wave.  The  time  dependent  phase  is 
C(i)  =  {ko  +  k)z{t)  -  (Of.  Electrons  separated  by  A^  =  2n  are  in  identical  positions 
in  adjacent  optical  wavelengths.  The  phase  velocity  is  y^o  =  L[{  kQ  +  k)%  -  k], 
and  the  dimensionless  time  is  x=ct/L.  The  energy  evolution  may  now  be 
written  as  y=  (e/cE/Ymc)cos(Co  +  vq  x).  This  can  also  be  restated  in  a 
dimensionless  time  regime  where  (...)  represents  a  dimensionless  time 
derivative  d{...)/clz,  so  that  y  =  (eEKL/ymc^)cos{(^  +  vqx).  Integration  of  this  term 
will  result  in 

Y  =  Yo  +  -^^^[sin(Co  +  Vol)  -  sin(Co)]. 
Yo  mc^o 

where  yo  is  the  initial  energy.  Rearranging,  the  time  varying  energy  in  the  weak 
optical  field  regime  is 

« 

7(^)  =  Yo  +  :^[  sin(;o  +  Vqx)  -  sin(Co)  ]  . 
Vo 

where  the  term  A  =  eKEL /yo  mc^  determines  the  energy  exchange.  This  shows 
that  the  energy  exchange  of  an  FEL  is  controlled  by  the  initial  energy  yo  of  the 
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electrons,  the  magnitude  of  the  electric  field  E,  the  magnetic  field  strength 
K  o^  B,  ant  the  electron  phase  Co  =  27czo/X  within  each  optical  wavelength,  and 
the  electron  phase  velocity  vq.  For  large  vo(vo»7c)  the  change  in  electron  energy 
oscillates  rapidly  along  the  undulator  and  is  small.  For  vq  near  resonance 
{vo  =  0),  the  electron  energy  changes  are  linear  in  z,  y(t)  =  yo  + /Axcos(Co).  and 
are  periodic  in  the  electron  phase  Co-  f^o''  -7c/2<Co<7c/2,  the  electron  energy 
increases,  and  for  7i/2<Co<37c/2,  the  electron  energy  decreases.  For  the  whole 
beam  randomly  spread  in  Co.  the  average  energy  change  is  zero.  To  this  order 
in  the  field  E,  there  is  no  gain. 

D.  THE  PENDULUM  EQUATION 

The  FEL  phase  and  phase  velocity  relationship  may  be  calculated  by 
returning  to  the  Lorentz  energy  expression,  y  =  -{e/mc)  '^^  ■  i?^.  FEL  gain  is  the 
net  energy  exchange  to  the  optical  field  and  must  be  calculated  to  higher  order 
in  the  field  E.  Substitution  of  the  optical  mode  electric  field  into  3-4  results  in 
an  expression  of  energy  gain  [8] 

Y       Tmc 

To  track  the  electron  phase  changes  as  the  energy  changes,  the  Lorentz  factor 
is  rewritten  y"^  =  1 -p^  =  1  -  (pf  +  p|)  where  ^i=  K/y.  For  highly  relativistic 
electrons  y  >  1 ,  and  (3^  =  1  -  (1  +  K^)'f.  Phase  velocity  C  =  v  =  (/cq  +  /()i  -  co, 
and  the  second  derivative  of  phase,  the  rate  of  phase  velocity  change,  is 

C  =  V  =  (/Co  +  k)z  =  kc  I^^^U  =  Kc  l^^^l  -f^  cos(C  +  (D).         3-5 

Equation  3-5  has  the  form  of  the  Pendulum  equation.   Taking  advantage  of  the 
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fact  that  the  initial  carrier  frequency  (Oq  =  /cqC  =  co  (1  +  K^)i2'f,  3-5  can  be 
restated  as  '< 

C  =  — 2 cos(C +  (!)). 

The  dimensionless  time  x=  c\iL,  wherein  0<t<  1  as  the  electrons  travel  down 
the  axis  of  the  undulator,  may  be  inserted  into  these  expressions.  This  requires 
that  d{.)i6t  =  L/C  d{..)/dx  or  (..)  becomes  (..),  and 

oo 

d%.)/dt^  =  {L/cf  d^{..)/dx^  =  (••)•  This  now  allows  the  Pendulum  equation,  3-5 
be  rewntten: 

oo  o      f  / 1^  2eEKkn 

^{'i)  =  y=i-\  — ^cos(C  +  (t))=  lalcos(C  +  (t)).  3-6 

where  the  dimensionless  term  \a\  =  AKNeKEL/yomc^  is  the  strength  or 
magnitude  of  the  optical  field  for  a  circularly-polarized  undulator  [9].  la  I  shows 
that  higher  velocity  (larger  y)  electrons  would  require  a  stronger  optical  field  E  in 
order  to  achieve  the  same  electron  bunching.  When  lal-cTt  the  field  is 
considered  to  be  in  a  "weak  field"  regime,  and  the  electron  bunching  after  a 
single  pass  through  the  undulator  is  negligible.  When  lal  =  7c  the  bunching 
effect  becomes  noticeable.  While  \a\:*>K  is  the  "strong  field"  case,  where  the 
electrons  quickly  bunch  within  the  undulator  and  may  cause  the  onset  of 
"trapped  particle  instability."  If  the  optical  power  is  strong  enough  to  cause  one 
synchrotron  oscillation  for  each  pass,  the  electrons  become  trapped  in  deep 
potential  wells.  This  results  in  beam  current  oscillation  that  destabilizes  the 
carrier  wave,  causing  sideband  frequencies  to  grow.  The  coherent  wave  that 
had  been  established  by  the  PEL  interaction,  may  be  amplified  or  degraded  by 
the  oscillating  current  of  the  trapped  electrons.  [9] 
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E.  THE  WAVE  EQUATION 

The  FEL  oscillator  starts  with  spontaneous  emission  and  creates  a 
classical  light  wave  that  is  amplified  by  stimulated  emission.  Because  the 
spontaneous  emission  bandwidth  of  the  wave  is  approximately  1/A/,  where  N  is 
the  number  of  periods  in  the  undulator,  there  is  some  measure  of  coherence  in 
the  initial  light  wave.  Since  A/»1 ,  the  bandwidth  is  sufficiently  narrow  to  form  a 
classical  wave  with  a  complex  field  envelope  {E{z,t),^{z,t))  whose  circularly- 
polarized  electric  and  magnetic  fields  are  expressed  by 

^ri^.t)  =  E(z.O[cos(v).-sin(y).0].     ^,(5?.0  =  E(z. f)[sin(v).cos(y).0]  3-7 

where  E{z,t)  is  the  wave  amplitude,  y  =  /cz  -  cof  -i-  ^{z,t),  and  the  phase  is  ^{z,t). 
The  amplitude  and  phase  vary  slowly  in  time  over  the  optical  period 
(E<t:coE,(i)<c(o<j)),  and  also  in  space  over  the  length  of  an  optical  wavelength  X, 
such  that  {E<s:kE,it:k(^)  where  k  =  2kJX  [8].  This  approximation  allows  the  fields 
to  be  derived  from  the  vector  potential  [9]: 

X{5?,0  =  -^^|^[sin(v),cos(y),0]. 


The  optical  wave  equation  becomes 


rfX(3?.f)  =  2 


aE  ^laE 

dz       c  dt 


[cos(\|/) -sin(v),0]  +  2E 


-4k 


dz       c  dt 


[-sin(\|/).-cos(\i/).0] 


:fi{^,t),  3-8 

where  uf  =  {V^  -  {Mc^)  {d^/dt^))  is  the  D'Alembertian  operator,  and  ^j^  is  the 

transverse  beam  current.     If  there   is  no  7^^  (no  transverse   motion  of  the 
electrons),  then  the  solution  to  3-7  is  E  =  E{z-ct)  and  ^  =  ^{z-ct). 
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Each  electron  contributes  to  the  total  transverse  cun-ent  so  that 

:7^  =  -ecZ^^8(3)(^-7^).  •  3-9 

where  7-  is  the  trajectory  of  the  /"'  electron  and  5^^'  is  a  3  dimensional  Dirac 
Delta  Function  [8].  The  coherence  volume,  V^c  is  the  product  of  the  gain 
bandwidth  length,  NX,  and  the  transverse  mode  area,  XL  =  XNXq,  so  that 
Vc  =  N^X^Xq.  In  order  to  correct  for  the  fast  variation  of  the  trigonometric 
functions  of  y  in  3-7,  an  average  is  taken  over  a  small  volume  element  dVc 
which  is  small  relative  to  the  coherence  volume,  but  large  relative  to  an  optical 
wavelength  X.  [9]   Using  this  averaging  and  inserting  3-7,  3-8  becomes 


dz      c  dt 


Ee''^  =  -AneKpiz  -  p^cf)<e-'^>  3-10 


where  p(z  -^z^t)  is  the  mean  density  of  the  electron  bunch.  If  the  longitudinal 
spread  of  the  bunch  is  so  great  that  the  density  is  no  longer  a  function  of  axial 
position,  there  will  be  uniform  evolution  of  the  electron  pulse  along  its  length, 
and  no  optical  modes  are  followed,  so  that  3-10  is  no  longer  dependent  on  z.  If 
the  number  of  undulator  periods  N  is  large,  there  is  a  small  bandwidth  and  y 
changes  only  slightly,  and  the  wave  equation  becomes 

a=-j<e-'^>  3-11 

where 

^  ^  AnNeKELe'^        .  ^  8N{enKLf)p 

The  field  magnitude  is  a,  and  j  is  the  dimensionless  beam  current  [8,9].  The 
beam  current,  y,  determines  the  reaction  of  the  optical  mode  to  the  bunching  of 
the  electron  beam,  and  provides  coupling  between  the  wave  and  the  electron 
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beam.  When  j  <  n  the  gain  is  low,  and  when  /  >7c  the  gain  is  high. 

For  a  linearly-polarized  undulator,  the  non-uniformity  of  electron 
acceleration  requires  that  these  expressions  be  modified  by  K^K{Jo(£,)  -  J^('^)). 
Here  K  is  the  undulator  parameter,  ^  =  K'^/2(1  +  K^),  and  Jq  and  J]  are  Bessel 
functions  of  the  first  and  second  kind.  [9] 

F.  DIMENSIONLESS  PARAMETERS 

Certain  dimensionless  parameters  recur  throughout  the  descriptions  of  the 
physics  of  the  FEL.  Their  use  allows  quick  summarization  of  performance 
without  regression  to  involved  calculations.  Several  important  dimensionless 
parameters  have  already  been  defined  and  are  used  in  this  paper.  They 
include  the  undulator  parameter,  K,  the  electron  phase  and  phase  velocity 
C  and  V,  the  beam  current  y,  and  the  optical  field  magnitude  la  I.  In  the 
derivation  of  3-10,  it  has  been  assumed  that  electron  and  optical  mode  sizes 
are  the  same.  The  filling  factor  F  =  nr^/nW^  is  the  ratio  of  the  area  of  the 
electron  beam  to  that  of  the  optical  mode  waist.  When  the  electron  beam  is 
entirely  within  the  optical  mode,  the  dimensionless  beam  current  j  (3-10)  should 
also  include  the  filling  factor  F  as  a  correction  for  the  smaller  volume  of  the 
electron  beam. 

The  transverse  dimensions  can  be  normalized  by  the  optical  mode  radius 


yn/LX.  Macroscopic  longitudinal  dimensions  can  be  normalized  by  the  undulator 
length  L.  The  normalized  electron  pulse  length  is  c52  =  IqINX  where  Ig  is  the 
electron  bunch  length.  The  Rayleigh  length  Zq  is  a  measure  of  optical 
diffraction  properties  and  is  dependent  on  optical  wavelength,  mirror  design  and 
configuration.   The  normalized  Rayleigh  length  is  zq  =  ZqIL.   Zq  is  related  to  the 
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optical  mode  waist  radius,  Wq,  by  Zo  =  iiW^/X.    The  normalized  mode  waist 


radius,    wq  =  Wo^nL/K    is    related    to    the    normalized    Rayleigh    length    by 


Wq  =  Wq<kLiX  =  ^XZqIk  <tiLK  =  -III. 

The  slippage  distance  As  is  the  difference  in  distances  traveled  by  the 
optical  wave  and  the  electron  beam  as  they  move  down  the  undulator.  This  is 
also  the  length  that  the  light  moving  at  c  passes  over  the  slower  electrons 
moving  at  p^c  during  this  period.  Their  relative  speed,  Av  =  {^  -  p^)c,  multiplied 
by  the  period  of  one  undulator  pass  At  =  Lie,  is  the  slippage  distance 
As  =  Av  At  =  {^  -  ^2)<^i^ic)  =  (1  -  Pz)^-  Inserting  the  definition  of  p^,  undulator 
length  L  =  NXq,  and  wavelength  X  =  XqO  +  K^)/2-f,  the  expression  for  slippage 
distance  becomes  As  =  NX.  The  normalized  electron  pulse  length  o^  =  Ig/NX  is 
scaled  by  the  slippage  distance  [9]. 
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IV.  THE  CEBAF  INFRARED  PEL  DESIGN 

A.   BACKGROUND 

The  Continuous  Electron  Beam  Accelerator  Facility  (CEBAF)  in  Newport 
News  Va.  is  administered  by  the  Southeastern  Universities  Research 
Association  (SURA)  for  application  to  nuclear  science,  technology  and  industry. 
The  CEBAF  superconducting  accelerator  incorporates  two  400  MeV  linear 
accelerators  (linacs)  with  interconnecting,  recirculating  electron  beam  lines. 
These  will  allow  five  simultaneous  electron  beams  at  energies  between  800 
MeV  and  4  GeV  to  be  delivered  to  three  end-stations.  A  program  is  underway 
to  design  and  construct  two  high  power  FELs  operating  at  the  near  infrared  (IR) 
and  the  deep  ultraviolate  (UV)  wavelengths.  The  design  configuration  of  the 
FEL  facility  will  allow  laser  operation  without  interference  to  the  nuclear  science 
research.  The  IR  FEL  will  take  advantage  of  the  operational  45  MeV 
accelerator  injector  as  its  driver,  while  the  UV  FEL  will  use  the  400  MeV  North 
Linac.  [10] 

The  CEBAF  facility  has  several  features  which  make  it  interesting  for 
FELs.  The  CEBAF  linacs  will  supply  the  FELs  with  a  high-quality  electron 
beam  with  normalized  emittance  £„  =  15  mm-mrad,  a  high  average  power  of 
800  kW,  energy  spread  ^yly=  .001,  and  a  duty  cycle  of  100  %.  The  fraction  of 
energy  extracted  by  the  FEL  is  1/A/,  so  that  the  average  power  of  the  FEL 
would  be  =  30  kW. 
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Figure  4-1 :  CEBAF  IR  and  UV  PEL  configuration.  The  UV  PEL  will 
utilize  the  400  MeV  North  Linac  for  its  driver,  while  the  IR  PEL  will  use 
the  45  MeV  Injector  [10]. 


These  are  important  parameters  because  if  the  PEL  were  ever  to  be  developed 
for  military  use,  high  average  power  and  performance  goals  would  have  to  be 
achieved,  and  this  experiment  will  be  an  important  step.  Additionally,  the 
accelerator  uses  advanced  superconducting  technology  that  will  be  incorporated 
into  the  next  generation  of  all-electric  drive  naval  vessels  that  are  now  being 
planned. 

The  Infrared  (IR)  PEL  is  a  multi-pass,  high-power  laser  operating  in  the 
near  infrared  at  a  nominal  wavelength  of  X  =  iSp-m  and  micropulse  frequency  of 
7.5  MHz  [1].  The  peak  current  is  72  A  with  a  short  micropulse  length  of 
/e  =  0.05  cm.    The   electron   beam   energy  is  {y-^)mc^  =  45   MeV  where  the 
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Lorentz  factor  is  y=89,  the  electron  mass  is  m,  and  c  is  the  speed  of  light. 
The  high-quality  electron  beam  has  energy  spread  of  only  Ay/y=  0.001  with  a 
normalized  emittance  of  e^  =  ye  =  15  mm-mrad.   [1 1] 

An  FEL  with  a  high-quality  electron  beam  operating  at  IR  wavelengths 
tends  to  have  an  electron  beam  that  is  much  smaller  than  the  optical  mode  in 
the  transverse  dimension,  and  smaller  than  the  slippage  distance  in  the 
longitudinal  dimension.  The  optical  mode  waist  radius  is  proportional  to  X^'^, 
and  the  slippage  distance  is  proportional  to  X.  In  the  transverse  direction,  mode 
distortion  tends  to  increase  the  gain,  while  in  the  longitudinal  direction,  the 
short-pulse  slippage  effects  tend  to  decrease  the  gain.  [11] 

In  the  CEBAF  case,  the  electron  beam  radius  is  r^  =  0.03  cm,  and  is  much 
smaller  than  the  radius  of  the  optical  mode,  lVo  =  0.15  cm.  Simple  theory 
indicates  gain  is  proportional  to  the  product  of  filling  factor  F  and  dimensionless 
beam  current  j,  but  the  actual  gain  observed  in  simulations  is  found  to  be  a 
function  of  filling  factor  F.  [10,1 1] 

B.   CEBAF  IR  FEL  PARAMETERS 

The  length  of  the  proposed  linearly-polarized  undulator  is  L  =  1.5  m  with 
wavelength  of  ?io  =  6  cm  over  N  =  25  periods.  The  linearly-polarized  undulator 
peak  field  strength  is  B  =  4.4  kG,  giving  an  undulator  parameter  of 
K  =  eBXoiZ^nmc^  ~  M^.  A  non-focussed  electron  beam  with  smooth 
propagation  and  minimum  electron  phase  velocity  spread  is  desirable  for 
injection  into  the  undulator  field.  This  requires  a  matching  condition  between 
the  rms  radial  beam  spread,  r,  and  the  rms  angular  spread  6,  such  that 
KRqT  =  V2Ye  [9].    The  emittance  is  e  =  r5,  so  that  the  matched  electron  beam 
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radius  is  then  r^  =  {yeko/^KKy^  =  0.03  cm.  The  off-axis  injection  caused  by  the 
radial  and  axial  spreads,  and  the  off-axis  focusing  fields  cause  a  slow, 
transverse  "betatron"  oscillation  of  the  electrons,  with  NK/y  =  0.5  betatron 
oscillations  along  the  undulator  [9].  The  amplitude  of  the  electron's  transverse 
wiggling  motion  is  Ax  =  KXq/^tc{  =  0.03  cm,  and  is  comparable  to  the  size  of  the 
electron  beam.  The  nominal  resonant  wavelength  is  X  =  Xo(1+K'^)/2/ =  15  ^im, 
but  the  CEBAF  IR  FEL  will  be  tunable  from  5M.m  to  20M.m.  The  dimensionless 
electron  beam  current  j  =  BN[KeKUJfpi-^mc^  ~  ^3^,  where  p  =  4x10^^  cm ~^  is 
the  electron  particle  density,  and  the  Bessel  function  factors, 
JJ  =  Jo(^Wi(^)  =  0.78  with  ^  =  K^/{20+K^))  =  0.38,  describe  the  reduced  coupling 
in  a  linearly-polarized  undulator.  [11] 

The  value  of  the  Rayleigh  length  that  minimizes  the  optical  mode  volume 
is  Zo=  /./VT2  =  43  cm  [9].  Assuming  the  optical  mode  is  shaped  as  a  lowest- 
order  Gaussian  with  its  waist  at  the  center  of  the  undulator  {xy,  =  0.5),  the  filling 
factor  averaged  along  the  undulator's  length  is  F  =  r^/[W§{^+L^/^2ZS)]~  0.027 
where  the  area  of  the  optical  mode  waist  is  kW§  =  ZqX  [9].  The  theoretical 
single-mode  gain  is  given  by  Gq  ~  0.135 jF  ~  50  %  [12],  and  does  not  include  the 
effects  of  beam  quality,  self-consistent  field  growth,  short-pulses,  or  optical 
diffraction.  When  the  interaction  includes  self-consistent  growth  of  the  optical 
field,  the  gain  increases  to  G  =  57  %.  The  electron  beam  energy  spread  has 
the  minor  effect  of  reducing  the  gain  by  only  0.5  %,  and  the  effect  of  emittance 
is  even  smaller.  The  CEBAF  IR  FEL  electron  beam  is  of  such  high  quality,  we 
will  not  include  the  effects  in  further  simulations.  [11] 
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C.  SHORT-PULSE  EFFECTS 

When  the  length  of  the  electron  pulse,  Ig,  is  comparable  to  the  slippage 
distance,  As  =  NX,  the  Fourier  components  of  the  pulse  current  are  comparable 
to  the  gain  bandwidth,  and  short-pulse  effects  influence  the  interactions  in  the 
undulator  [9].  As  an  optical  pulse  bounces  off  the  resonator  mirrors  and  arrives 
at  the  entrance  to  the  undulator,  at  i  =  0,  an  electron  micropulse  from  the 
accelerator  is  timed  to  arrive  simultaneously.  A  slight  displacement  of  one  of 
the  resonator  mirrors,  termed  the  "desynchronism"  d  =  As/nX,  alters  the 
distance  the  rebounding  light  must  travel,  and  therefore  will  serve  as  a  means 
to  adjust  the  arrival  of  the  optical  pulse  relative  to  the  subsequent  electron 
pulses.  Desynchronism,  d,  is  normalized  to  the  slippage  distance,  and  an  exact 
synchronism  between  arrivals  of  the  light  and  the  electron  pulses  occurs  when 
d  =  0. 

For  the  CEBAF  IR  FEL  operating  at  the  nominal  wavelength  X=^5\l^^, 
Figure  4-2  shows  an  undulator  with  exact  synchronism  d=0.  The  window  for 
calculation,  at  middle  left,  is  four  slippage  distances  long,  and  travels  with  the 
light  pulse  at  speed  c.  The  electron  pulse,  traveling  at  slower  speed  p^c,  slips 
backward  through  one-fourth  of  the  window  per  pass.  The  positions  of  the 
short  electron  pulse  relative  to  the  optical  pulse  is  shown  at  left  bottom.  The 
electron  pulse  position  y(z-x),  initially  is  shown  as  the  lighter  parabolic  shape  at 
T  =  0.  As  bunching  develops  the  electron  pulse  slips  backward,  and  the  final 
position  is  indicated  by  the  darker  shape  at  t=  1.  The  pulse  form  is  parabolic 
shape  j{z)  =  y(1-2z^/o|)  for  y(z)>0,  or  otherwise  zero.  The  electron  normalized 
pulse  length  o^  =  lg/NX=  1.4,  and  the  peak  current  is  jF  =  3.8  where  F  =  0.027. 
The  optical  amplitude  grows  proportionally  to  j  on  each  pass,  while  amplitude 
losses  due  to  the  resonator  are  calculated  by  e"^'°  for  each  pass.  [9,1 1] 
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Figure  4-2:  Short-pulsa  effects  in  an  undulator  with  exact  synchronism 
{d  =  Q). 
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In  order  to  start  the  optical  pulse  growth,  random  shot  noise  Is  generated 
by  addition  of  a  small  phase,  5C  =  0.0001 ,  to  each  of  the  sample  electron 
phases.  The  optical  pulse  is  free  to  evolve  subject  to  gain  effects  of  the 
electron  pulse,  and  resonator  losses  due  to  Q  over  n  =  800  passes  through  the 
undulator.  Each  pass,  a  new  electron  enters  the  undulator  to  drive  the 
rebounding  optical  pulse.  The  shot  noise  creates  a  low  power,  broadband 
optical  pulse  in  the  first  few  passes  with  a  length  that  is  roughly  equal  to  the 
sum  of  the  slippage  distance  and  the  electron  pulse  length,  o^  +  1.  This  makes 
sense  because  this  is  the  end-to-end  distance  over  which  the  optical  pulse  is 
exposed  to  the  amplifying  electrons  when  at  perfect  synchronism. 

The  optical  field  is  amplified  at  the  rear  of  the  optical  pulse  centroid  after 
n  =  10  passes.  The  electrons  are  amplifying  the  optical  pulse,  which  has 
partially  passed  over  the  electrons.  It  is  clear  that  as  the  number  of  passes,  n, 
increases,  the  centroid  of  the  optical  pulse  moves  backward  (to  the  left)  due  to 
the  continual  amplification  at  the  rear  of  the  pulse.  The  amplitude  evolution 
peaks  near  n  ~  200  passes,  and  then  decreases.  The  final  shape  of  the  pulse, 
\a{z,n)\,  is  shown,  top-left,  above  the  amplitude  evolution  plot.  The  maximum 
amplitude  of  the  optical  field  is  \a{z,n)\  =  58  is  shown  on  the  grey  scale  in 
black.  A  single  contour  line  signifies  half-maximum  values,  while  zero 
amplitude  is  in  white. 

The  power  spectmm  evolution  (middle-center)  starts  near  peak  gain  at 
vq  =  2.6,  and  grows  in  width  at  n  =  150  passes.  The  broadened  range  of  phase 
velocities  v  at  this  point  is  caused  by  a  short  pulse.  The  weak-field  gain 
spectnjm,  G(v),  is  plotted  on  the  same  scale  below  for  reference.  The  final 
power  spectrum,  P(v,n),  shows  the  spread  in  phase  velocities  for  the  last  pass, 
n  =  800.   The  square  tick-mark  at  the  top  of  the  power  spectrum  indicates  the 
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center  of  the  final  power  spectnjm,  while  the  pointed  tick-mark  shows  the 
central  wavelength  of  the  initial  radiation  at  resonance.  The  power  evolution, 
P{n),  also  shows  rapid  increase,  peaking  near  n  =  150,  then  tapers  ultimately  to 
zero  power.  The  light  wave  in  exact  synchronism  systems  eventually  drifts 
away  from  the  slower  electrons  leaving  the  FEL  with  steady-state  power  of 
zero. 

The  electron  beam  phase  velocity  evolution  (right-center)  shows  a  single 
phase  velocity  emerging  because  of  the  vanishing  optical  power.  On  the  final 
phase  velocity  plot,  f{v,n)  (right-top),  the  coincident  initial  and  final  phase 
velocity  tick-marks  indicate  near-constant  energy  of  the  electron  pulse  due  to  a 
optical  field. 

In  an  undulator  with  desynchronism  d  >0,  the  centroid  of  the  optical  pulse 
arrives  at  x=0  slightly  ahead  of  the  centroid  of  the  electron  pulse  which  then 
passes  back  over  the  optical  pulse.  As  the  pulses  propagate  and  x  increases, 
the  electrons  bunch  and  radiate,  amplifying  the  rear  portions  of  the  optical  pulse 
with  gain.  The  shape  of  the  optical  pulse  is  distorted  with  each  successive 
pass  through  the  undulator.  Because  the  light  pulse  essentially  grows  from  the 
rear,  the  light  pulse  will  be  distorted  in  amplitude  over  the  slippage  distance  as 
the  rear  of  the  pulse  grows.  Figure  4-3  shows  this  effect  with  a  small 
desynchronism  d  =  0.005,  for  the  CEBAF  IR  FEL  described  in  Fig  4-2. 

In  Fig  4-3,  altering  the  FEL  with  small  desynchronism  d  =  0.005  results  in 
major  changes  is  performance  of  the  laser  system.  The  optical  amplitude 
evolution  survives  and  produces  an  amplitude  in  steady  state  at  n  =  300  passes 
which  has  the  dual-peak  modulation  in  the  final  pass  window  \a{z,n)\. 
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Figure    4-3:    An    FEL    system    with    desynchronism    and    amplitude 
modulation  due  to  short-pulse  effects. 
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The  power  spectrum  evolution,  P{v,n),  and  optical  power,  P{n),  also 
achieve  steady  state,  with  the  maximum  power  much  greater  for  this  laser  with 
non-zero  desynchronism  than  the  earlier  synchronous  FEL. 

Though  the  final  optical  pulse  shape  \a{z,n)\  is  frozen  in  a  single  picture, 
the  series  of  pictures  for  each  successive  pass  shows  the  broad  amplitude 
subpulses  originate  in  an  area  of  z  =-1,  and  grow  while  moving  to  the  right. 
This  effectively  describes  the  optical  pulse  being  amplified  along  its  rear,  and  a 
subpulse  passing  forward  away  from  the  slower  electrons.  The  magnitudes  of 
the  subpulses  would  depend  on  the  parameters  of  the  system  under  simulation. 
The  amplification  and  growth  of  the  rear  of  the  optical  pulse  also  has  the  result 
that  the  centroid  of  the  pulse  moves  at  a  speed  less  than  c,  which  effectively 
lowers  the  pulse  group  velocity. 

The  tick  marks  in  the  final  power  spectrum  picture  shows  the  spread  of 
center  phase  velocities  expected,  and  there  is  a  small  sideband  caused  by 
trapped-particle  instability.  The  wider  spread  of  the  electron  spectrum  f(v,n)  is 
proportional  to  the  separatrix  peak-to-peak  height,  and  the  initial  and  final  tick 
marks  show  a  drop  in  electron  phase  velocities,  indicating  a  net  transfer  of 
energy  to  the  optical  field. 

The  optical  wavelength  of  an  FEL  is  defined  k  =  Xo((1  +  K^)/2-f).  In  order 
to  tune  an  FEL  to  a  new  wavelength,  the  undulator  parameter,  K  =  eBXQ/2nmc^, 
is  altered  by  changing  the  mean  magnetic  field  E.  This  in  turn  changes  the 
optical  mode  waist,  Wq,  the  filling  factor,  F,  the  dimensionless  electron  beam 
current,  y,  the  dimensionless  pulse  length,  o^  =  Ig/NX,  and  several  other  FEL 
parameters.  The  effects  short  wavelengths  for  the  CEBAF  IR  FEL  operating  at 
X  =  5|im  is  shown  in  Fig  4-4. 
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Figure  4-4:  An  FEL  operating  at  X  =  5M.m  with   long  electron  pulse 
lenqth  a.. 
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The  dimensionless  pulse  length  here  is  o^  =  4.3,  three  times  the  length  of 
the  electron  pulse  for  the  FEL  operating  at  A.=  ISjim  shown  in  Fig  4-3.  The 
amplitude  evolution  plot  clearly  shows  amplitude  spikes  growing  from  the  rear  of 
the  optical  pulse,  and  moving  forward  over  several  passes.  The  longer  pulse 
results  in  a  longer  optical  pulse  spread  along  eight  or  ten  slippage  distances 
and  \a{z.n)\  amplitude  spikes  spaced  roughly  one  slippage  distance  apart.  In 
contrast  to  the  narrow  power  spectnjm  of  the  short  pulse  FEL  of  Fig  4-3,  here 
the  power  spectrum  is  spread  among  several  sidebands  without  a  steady-state. 

Power  evolution,  P{n),  rises  quickly  but  never  achieves  a  steady  value. 
The  higher  average  power  is  a  direct  result  of  the  longer  electron  pulse 
interacting  with  the  optical  pulse  for  a  longer  period  and  so  transferring  more 
energy.  The  electron  spectrum,  /(v.n),  of  Fig  4-4  is  again  broad,  with  the 
average  phase  velocity  dropping  due  to  net  energy  transfer  from  the  electrons 
to  the  optical  pulse. 

While  electron  pulse  described  by  o^  =  4.3,  shown  in  Fig  4-4  has  the 
smallest  slippage  within  the  tuning  range  5)im  ^  X  <  20M,m,  the  largest  slippage 
is  o^  =  1.1  at  20M,m.  This  system  is  shown  in  Fig  4-5.  Figures  4-4  and  4-5 
serve  well  to  contrast  the  various  effects  of  pulse  length  on  system 
performance.  Short-pulse  effects  generally  result  in  short,  narrow  optical  pulses 
spread  over  fewer  slippage  distances  and  comprised  of  fewer  spikes. 

Figure  4-6  shows  an  FEL  operating  near  X=15)im  with  relatively  high 
resonator  losses,  0  =  50,  and  a  moderate  desynchronism  d  =  0.0^65.  This 
higher  value  of  d  causes  the  modulation  of  the  amplitude  evolution  to  stand  out 
cleariy.  [11] 
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Figure  4-5:  FEL  operating  aXX  =  20p.m  and  short  electron  pulse  length 


a^  =  1.1, 
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Figure  4-6:  An  FEL  with  steady  state  limit  cycle  behavior  in  the  power 
evolution  P{n).  [1 1] 
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The  power  evolution  P{n)  exhibits  sharp  overshoot  followed  by  a  steady 
state  oscillation  in  power  as  new  optical  pulses  grow  periodically.  The  period  of 
the  oscillations  appears  to  be  An  =  90  passes.  This  oscillation  is  known  as 
"limit-cycle  behavior."  This  phenomena  occurs  in  regimes  of  moderate 
desynchronism  where  a  stable  overlap  does  not  occur  between  the  optical 
pulse  and  the  electron  pulse.  Moderate  field  strengths  cause  about  one 
synchrotron  oscillation,  which  causes  the  optical  pulse  to  continually  shift  from 
one  shape  to  another.  [9,1 1] 

D.  CEBAF  IR  FEL  TUNING  RANGE  AND  GAIN 

The  tuning  range  of  the  CEBAF  IR  FEL  (5)im  <,  X^20M.m),  combined  with 
the  effects  of  varying  desynchronism  on  gain,  allows  for  a  broad  range  of 
system  performance.  The  broad  tuning  range  of  the  FEL  includes  electron 
pulse  lengths  ranging  ^.^<C2  <4.3,  which  would  imply  a  wide  range  of  system 
gain.  Neglecting  the  effects  of  the  slight  energy  spread,  and  using  multiple 
passes  to  achieve  steady-state  gain,  simulations  were  mn  over  a  range  of 
desynchronisms. 

The  shortest  wavelength  of  the  tuning  range,  X  =  5p.m,  has  the  longest 
pulse  length  o^  =  4.3.  Varying  desynchronism  d  from  0.005  through  0.06 
resulted  in  a  steady  single-pass  gain  of  G  =  170  %.  For  wavelength  X=  10)j.m, 
with  pulse  length  o^  =  2.2,  the  gain  is  only  about  60  %.  At  X  =  15|j.m,  the  gain  is 
about  35  %.  At  the  longest  wavelength  of  this  tuning  range,  X  =  20|im,  the  gain 
is  about  20  %.  Fig  4-7  shows  the  variation  of  the  system  gain  across  the 
wavelength  tuning  band  of  the  CEBAF  IR  FEL. 
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Figure  4-7:  Gain  as  a  function  of  wavelength  X  for  the  CEBAF  IRFEL 


The  longest  wavelength  X  =  20|im  has  the  lowest  gain  in  the  tuning 
bandwidth  with  G  =  22  %.  We  considered  lengthening  the  undulator  from 
A/  =  25  to  A/  =  35  periods  in  order  to  increase  gain.  The  dimensionless  current  / 
increased  from  138  to  379.  The  undulator  length  L  increased  from  150  cm  to 
210  cm,  while  the  dimensionless  electron  pulse  length  decreased  from  a^  =  1.1 
to  0.77.  The  average  filling  factor  F  was  reduced  by  half,  but  since  ;  increased, 
the  product  jF  increased  by  a  factor  of  three.  While  increasing  jF  increases 
gain,  the  decrease  in  a^  decreases  gain.  The  net  result  from  the  lengthening  of 
the  undulator  is  determined  by  simulation.  The  gain  is  improved  to  G  =  36  % 
for  /V  =  35  periods  and  G  =  50  %  for  /V  =  45  periods.  The  cost  of  increasing  the 
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length  of  the  undulator  in  order  to  improve  gain  at  longer  wavelengths  would  not 
be  large. 
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V.  OPTICAL  MODE  DISTORTION 

A.  BACKGROUND 

After  originating  from  spontaneous  emission,  the  optical  field  in  the 
undulator  is  a  sum  of  rebounding  light  waves  superimposed  with  newly 
generated  light  from  the  bunched  electrons  within  the  electron  beam  [15].  The 
longitudinal  effects  of  short  pulses  were  described  in  Chapter  IV,  and  the 
transverse  effects  caused  by  natural  diffraction  and  distortion  by  the  electron 
beam  will  be  treated  in  this  chapter.  The  effects  of  short  pulses  will  not  be 
included  in  this  chapter.  For  the  CEBAF  IR  PEL  T  =  0.027  [10].  As  electron 
bunching  develops  late  in  the  length  of  the  undulator,  the  new  light  is  created  in 
an  area  that  is  =  F  times  the  size  of  the  mode  area.  The  newly  created  light 
from  the  small  electron  beam,  has  a  short  Rayleigh  length  and  diffracts  rapidly. 

B.  GENERAL  MODE  DISTORTION 

The  effects  of  the  small  electron  beam  size  acting  on  the  optical  mode 
may  be  simulated  through  numerical  integration  of  the  parabolic  wave  equation, 
3-9,  coupled  with  the  Lorentz  force  equation,  3-1.  To  calculate  the  effects  of 
mode  distortion,   it  is   necessary  to  shift  the  transverse  dimensions  into  a 


dimensionless  coordinate  system,  normalized  to  the  mode  size,  x  =  x^ki2L . 
The  optical  amplitude,  la  I,  and  the  dimensionless  current  density,  y,  are  also 
shifted  into  this  coordinate  system.  The  parabolic  wave  equation  now  becomes 


_1^2  .    a 


^dx 


a{it,x)  =  ^ie''^>,  5-1 
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where  the  operator  ^^  =  d^/dx^  +  d^/dy^,  properly  describes  the  diffraction  of  the 
optical  wave  in  transverse  dimensions  [13].  An  important  effect  of  diffraction  is 
self-focusing  of  the  optical  mode  back  into  the  electron  beam.  The  optical 
phase  shift,  ^{x),  occurs  in  the  center  of  the  optical  mode,  near  the  electron 
beam,  and  acts  to  focus  the  light  back  onto  the  undulator  axis.  The  FEL 
interaction  causes  the  light  to  remain  closer  to  the  electron  beam  than  would 
occur  in  the  presence  of  natural  diffraction  alone.   This  effect  is  greater  when 


gain  is  larger.    When  the  optical  mode  is  larger  than  szUk    the  distortion 
operator  ~  1,  and  can  be  neglected.  [13] 

C.  MODE  DISTORTION  SIMULATIONS 

Figure  5-1  shows  a  numerical  simulation  of  a  Gaussian  optical  wavefront 
propagating  without  the  effects  of  electron  beam  current  (y  =  0),  as  x  goes  from 
0  to  1  in  an  undulator  with  N  =  25  periods.  The  wavefront  has  weak  initial 
optical  field  amplitude,  ao=l.O,  at  mode  center,  and  mode  curvature  giving 
normalized    Rayleigh    length    of   Zq  =  0.5    focussed    at   x^  =  0.5.     Transverse 


dimensions  are  normalized  to  the  mode  waist  <ki2L  .  The  grey  scale  at  upper 
right  shows  the  intensity  for  minimum  (white),  and  for  maximum  (black),  with 
five  white  contour  lines  to  illustrate  the  focusing  and  diffraction  of  the  light  wave. 

The  optical  amplitude  evolution  plot,  la(x,x)l,  (upper-left)  shows  the  optical 
mode  evolution  over  a  single  pass  through  the  undulator.  The  five  contour  lines 
show  the  focusing  and  diffraction  expected  for  an  undistorted  optical  wavefront. 
The  optical  amplitude  is  greatest  at  the  mode  center,  and  decreases  radially  as 
shown  in  the  transverse  optical  amplitude  plot,  la(x,y)l,  representing  the  final 
mode  at  the  end  of  undulator.  The  optical  phase  shift,  (|)(0,x),  shifts  negatively, 
A<t)  =  -x/Zq,  during  diffraction. 
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Figure  5-2  shows  an  FEL  with  the  same  optical  parameters  as  shown  in 
Fig  5-1,  but  with  optical  mode  distortion  caused  by  an  electron  beam  current 
y  =  138  for  the  CEBAF  IR  FEL.    The  electron  beam  has  normalized  radius, 


o^  =  rt,HLX/K  =  0.125,  and  has  initial  phase  velocity  vo  =  4.6.  With  a  single 
contour  line,  the  optical  amplitude  plot,  la(x,T)l,  shows  focusing  of  the  light  with 
Rayleigh  length  zq  =  0.5  in  toward  the  mode  waist  positioned  at  t^  =  0.5. 
Relative  to  the  undistorted  field  of  Fig  5-1,  the  field  distortion  in  Fig  5-2  is  clear. 
The  maximum  amplitude,  la(x,T)l,  is  greater,  and  the  contour  line  spreads 
rapidly  from  the  beginning  of  the  undulator  as  the  distortion  of  the  electron 
beam  drives  up  the  amplitude  of  the  optical  mode  along  the  second  half  of  the 
undulator.  The  final  transverse  amplitude  plot,  la(x,yl,  shows  the  relatively 
narrow  mode  focussed  around  the  electron  beam.  The  electron  bunching 
current,  o(x,t)  «  y(r)<cos(^,(t))>,  evolves  as  the  electron  phases,  C.  bunch  within 
each  optical  wavelength.  The  plot  of  final  electron  bunching,  c{x,y),  shows  the 
relative  disparity  in  size  of  the  electron  beam  relative  to  the  final  optical  mode, 
la(x.y)l.  [11] 

The  electron  phase  velocity  evolution,  /(v.t)  in  Fig  5-2,  starts  at  vq  =  4.6, 
and  broadens  slightly  due  to  bunching.  There  is  also  slight  bunching  shown  in 
the  phase-space  plot,  (C,v).  The  optical  phase  shift,  ^{0,x),  is  also  affected  by 
distortion  effects  of  the  beam  current,  j.  From  5-1 ,  the  FEL  optical  phase  shift 
0  =  y<sin(C  +  (t))>/lal,  is  maximized  at  (C  +  <t))  =  7c/2.  The  electron  phase,  C.  is  a 
function  of  dimensionless  time,  x,  so  that  for  strong-field  FEL  interaction, 
(|)  =  yi/la  I,  and  for  weak  fields,  ^  ~  yT^/12.   [9] 
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FEL   WAVEFRONTS 
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Figure  5-1 :  Non-distorted  optical  mode  propagation  diffraction. 


In  either  a  strong  or  weak  fields,  the  FEL  interaction  gives  a  positive  shift 
to  the  optical  phase.  Variation  of  4){0,t)  as  x  increases  from  0  to  1  is  driven  by 
electron  beam  interaction. 

Exponential  growth  of  gain  and  power  is  delayed  in  the  undulator  until  the 
bunching  time,  xg  =  {2/j)'^'^  =  0.25,  has  passed,  after  which  the  strong  beam 
current,  y,  drives  up  the  gain  and  power  [9].  The  gain  for  the  FEL  in  Fig  5-2  is 
ln(i  +  G{t))  =  0.846  or  G  =  133%.  Theoretical  single-mode  gain  for  this  system 
predicts  only  Gq  =  0.135;^  =  43%.  in  single-mode  theory,  the  dimensionless 
current  density,  J,  is  proportional  to  the  electron  beam  density,  p,  and  is 
therefore  inversely  proportional  to  the  square  of  the  electron  beam  radius,  or 


50 


FEL  WAVEFRONTS 


2   ia(x,t)  i 


2 

10 


a  (:'.,!) 


a{x,y)  : 


c^iy-r/} 


# 


j=138 

X  =0 

o 

V   =4.6 
o 

T   =0.5 
w 


a  =0.125 

X 

a  =1 
o 

z   =0.5 
o 

N=25 


0  |a(x,X)|  2 

0  |a(x,y)|  2 

0  a(x,T)  70 

0  a(x,y)  70 


f  (v.- ) 


-10 


wm 


1  -rt/2        I         371/2  0 


Figure  5-2:  FEL  with  optical  mode  distortion.  [11] 
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The  filling  factor,  F,  is  proportional  to  the  electron  beam  radius,  F  «  r/, 
and  so  the  theoretical  gain,  Gqoc  jF,  becomes  independent  of  the  electron 
beam  size  as  long  as  the  electron  beam  is  contained  within  the  optical  mode. 
But  single-mode  theory  does  not  account  for  the  mode  distortion  caused  by  the 
intense  current  and  small  filling  factor  of  the  FEL  in  Fig  5-2.  [1 1] 

Figure  5-3  shows  the  effects  of  shifting  mode  waist  position,  x^,  and 
Rayleigh  length,  zq,  to  achieve  increased  gain.  Showing  strong  diffraction  and 
distortion  effects,  the  amplitude  evolution,  la(x,T)l,  shows  the  light  focussed 
inward  to  x^  =  0.19,  and  then  sharply  diffracting  outward  from  a  narrow  waist 
near  x  =  0.40.  This  is  due  in  part  to  the  short  Rayleigh  length,  zq  =  0.14,  which 
causes  light  to  diffract  rapidly.  The  plot  of  the  final  mode,  la(x,y)l,  clearly 
shows  the  narrow  focusing  of  the  optical  mode  around  the  electron  beam,  and 
the  rapid  spread  of  highly  diffracted  light. 

The  larger,  irregular  shift  of  ^(O.t)  for  the  FEL  in  Fig  5-3  is  due  to  the  short 
Rayleigh  length,  zo  =  0.l4,  and  the  varying  effects  of  the  field  amplitude  la! 
along  the  undulator.  The  apparent  discontinuity  late  in  the  plot  of  ^(0,x) 
represents  the  optical  phase  shift  decreasing  to  =  -k.  Electron  bunching 
occurs  later  in  the  undulator  than  for  the  FEL  of  Fig  5-2,  as  shown  by  the 
bunching  current  plot,  o(x,t).  The  final  gain  is  much  greater  than  for  Fig  5-2, 
with  ln(1  +  G(t))  =  1.14,  or  G(i)  =  2l2%.  Single-mode  theory,  neglecting  the 
effects  of  mode  distortion,  predicts  this  FEL  to  have  gain  Go  =  40  %. 
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Figure  5-3:  An  FEL  with  short  Rayleigh  length  and  high  distortion. 
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D.  MODE  PARAMETERS,  RESONANCE  AND  GAIN 

Variation  of  Rayleigh  length  and  mode  waist  position  and  the  resultant 
gain  curves  are  shown  in  Figures  5-4  and  5-5.  Typical  values  for  Rayleigh 
length  and  mode  waist  position  are  zq  =  nW^/LX  =  Zq/L  =  I/VT2  or  zq  =  0.5,  and 
Xjv  =  0.5,  but  simulations  for  the  CEBAF  IR  PEL  produced  peak  gain  at 
zq  =  0.33.  Using  this  Rayleigh  length,  the  peak  gain  for  variation  of  mode  waist 
position  was  located  at  x^  =  0.35. 

Figure  5-6  shows  the  optical  gain  spectrum,  G(vo),  plotting  the  gain 
against  the  PEL  phase  velocity,  vo  = /.[{/c+/co)p2-/c]  [7].  Rayleigh  length, 
zo  =  0.14,  and  mode  waist  position,  Xy^  =  0.^9,  are  selected  to  maximize  the 
single-pass  gain  for  the  small  CEBAP  electron  beam.  When  X=15|im,  the 
CEBAP  IR  PEL  is  near  resonance  vq  =  0,  and  Pig  5-6  shows  the  gain  around 
resonance.  When  the  phase  velocity  is  off-resonance  by  Avq  =  n,  the  shift  in 
wavelength  is  aX/X  =  Avq/2kN  =1/2/V  =  2  %.  The  electron  beam  has  a  radial 
parabolic  shape  described  by  j{r)  =  j{^-^^/2og)  for  r  <  V2ae,  ^rid  zero  othenwise. 
The     peak    current     is    j  =  138,    and    the     normalized    beam     radius    is 


Og  =  rg^nJLX  =  0A25.  [11] 

The  simulation  in  Pig  5-6  calculates  gain  including  the  transverse  effects 
of  diffraction,  but  ignores  the  longitudinal  short-pulse  effects.  At  each  value  of 
vo,  a  wavefront  enters  the  PEL  undulator  with  peak  field  strength  ao  =  1  in  the 
weak  field  regime,  where  ao  =  AnNeKLEz-fmc^,  and  E  is  the  optical  electric  field 
strength.  The  peak  gain  in  Pig  5-6  is  G  =  233  %  at  phase  velocity  vq  =  5.  The 
region  of  positive  gain  is  larger  than  the  region  of  negative  gain,  as  is  typical  of 
high  current  PELs  [9].  The  gain  spectmm  tails  off  smoothly  as  vq  increases  on 
the  right  of  Pig  5-6.  [11] 
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Figure  5-4:  Variation  of  gain  with  Rayleigh  length,  Zq. 
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Figure  5-5:  Variation  of  gain  with  mode  waist  position,  z^. 
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***  FEL  3d  simulation,  single-pass  gain  *** 
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Figure  5-6:  The  optimized  CEBAF  IF  FEL  gain  spectrum  including  only 
transverse  modes. 
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E.  GAIN  AND  SMALL  ELECTRON  BEAMS 

The  FEL  gain  observed  for  Figs  5-2  and  5-3  indicates  that  single-mode 
theory  can  significantly  underestimate  the  effects  of  diffraction  and  distortion  on 
actual  FEL  gain.  Calculating  the  single-pass  gain  over  a  broad  range  of  filling 
factors,  the  dependence  of  FEL  gain  on  electron  beam  size  is  examined  in  Fig 
5-7.  Neglecting  optical  mode  distortion,  the  single-mode  FEL  gain  is 
Go  =  0.135yr.  The  gain  ratio  G/Gq  is  determined  numerically  and  plotted  as  a 
function  of  the  electron  beam  size  compared  to  the  optical  mode  waist,  Og/wo. 
with  a  single  line  on  the  plot  signifying  G/Go  =  l.O.  The  product  of  the  peak 
beam  current  and  the  beam  area,  ya| «  jF,  is  held  constant  at  jF  =  0.1.  The 
electron  beam  shape  is  taken  to  be  parabolic,  j{r)  =  y(1-r^/2o|)  for  r  <  V2oe  and 
zero  othenwise.  A  small  value  of  jF  is  used  to  insure  that  there  are  no  high 
gain  effects  for  values  on  the  left  side  of  Fig  5-7  where  F  =  Cg/w§  is  small  and  J 
may  be  large.  The  Rayleigh  length  is,  Zq  =  1/VT2  =  0.29,  for  optical  single-mode 
gain  Gq  and  position  of  the  mode  waist  is,  x^^  =  0.5.  These  values  are  chosen 
to  describe  a  more  typical  FEL  mode  instead  of  the  optimized  CEBAF  IR  FEL. 
The  phase  velocity  is  vq  =  5  for  peak  single-mode  gain  [12],  and  the  peak 
optical  amplitude  is  ao  =  1  for  weak  fields.  [11] 

The  plot  shows  a  clear  dependence  on  the  electron  beam  size  at  all 
values  plotted.  The  single-mode  theory  is  correct  for  the  range  of  values  Og/wo 
where  the  G  =  Gq.  We  see  that  there  is  no  range  of  values  where  the  single- 
mode  theory  is  correct,  but  only  a  single  point.  For  a  small  electron  beam,  the 
single-pass  gain  including  many  modes  can  be  significantly  higher  than  the 
single-mode  gain,  Gq.  It  is  emphasized  that  the  mode  distortion  crucial  to  the 
increased  value  of  G/Gq  is  not  caused  by  high  current.  [11] 
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***  FEL  3d  simulation,  single-pass  gain  *** 
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Figure  5-7;  General  gain  curve  for  FELs  with  small  electron  beams. 
[11] 
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When  the  product  jP  is  decreased  by  xioo  to  jF  =  0.001,  the  plot  in  Fig  5- 
7  remains  substantially  the  same.  When  jF  is  increased  to  unity,  or  above,  the 
gain  enhancement  due  to  mode  distortion  is  further  increased  by  high-gain 
effects  on  the  left  side  of  Fig  5-7  where  beam  radius  is  small  and  j  is  large.  It 
is  not  clear  from  this  analysis  what  part  of  the  distorted  mode  will  remain  in  a 
resonator  after  many  passes.  The  actual  gain  of  an  FEL  oscillator  may  be 
reduced  from  that  shown  in  Fig  5-7  depending  on  the  resonator  design 
[13,15,16].  No  known  work  explains  Fig  5-7,  and  in  particular,  its  independence 
of  the  value  of  jF  below  unity,  where  there  is  vanishingly  small  mode  distortion. 

[11] 

F.  PROPAGATION  AND  INTERFERENCE 

A  modification  of  the  simulation  that  produced  Figures  5-1  to  5-3  subtracts 
the  input  optical  field  from  the  total  field  to  allow  analysis  of  the  propagation  of 
newly  created  light  during  a  single  pass  through  the  undulator.  This  process  is 
interesting  because  of  the  complex  optical  field  interactions  shown  in  Fig  5-8. 
Because  the  amplitude  plots,  Ala(x,x)l  and  Ala(x,y)l,  are  generated  by 
subtraction,  it  is  possible  to  have  negative  amplitude  values,  representing 
destnjctive  interference  in  the  total  optical  field,  reducing  field  amplitude. 
Destructive  interference  occurs  when  "new"  light  from  the  bunching  of  electrons 
encounters  "old"  light  with  phase  difference  of  \A(^\  =  AxIzo  =  k  radians.  Noting 
that  the  newly  generated  light  appears  in  the  amplitude  evolution,  Ala(x,x)l,  near 
Ax  =1/3,  and  holding  optical  mode  waist,  wq  =  [2o+ 1/l2zo]^^,  constant  by  the 
input  Rayleigh  length,  zq  =  w^  =  Cg,  a  phase  shift  relationship  may  be  found: 
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OPTICAL  MODE  DISTORTION: PROPAGATED 
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Figure  5-8:  Negative  interference  of  propagated  light. 
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For  electron  beam  radius,  Og  <  0.3,  conditions  for  destnjctive  interference  may 
exist. 

The  plot  of  propagated  optical  amplitude  evolution,  Ala(x,T)l,  shows  an 
area  of  destructive  interference  starting  near  x  =  0.7,  which  is  obvious  in  the 
final  amplitude  plot,  Ala(x,y)l,  as  a  light  area  in  the  center  of  the  optical  mode. 
To  evaluate  the  interference  effect  over  a  range  of  values  of  yo|,  the  product 
was  held  constant  while  y  and  Cg  were  varied.  Destructive  interference  was 
observed  for  0.12  <Oe<  0.3,  and  the  point  of  initial  destructive  interference 
occurred  earlier  in  the  undulator  as  the  electron  beam,  Cg,  is  narrowed.  This  is 
in  part  caused  by  the  constant  product  jog,  which  rapidly  raises  beam  current, 
j,  as  the  electron  beam  narrows.  This  appears  as  a  radial  'dimple'  across  the 
propagated  wavefront  after  a  single  pass  through  the  undulator.  In  order  to 
determine  whether  the  phenomenon  of  destnjctive  interference  in  the  optical 
wavefront  continues  outside  the  undulator,  a  multi-pass  simulation  within  the 
entire  resonator  is  required.  This  destructive  interference  over  a  single  pass 
was  first  investigated  by  LT  Carl  Bice  in  his  thesis  in  December,  1991,  and  is  a 
continuation  of  that  work.[14] 


61 


VI.  MULTI-PASS  OPTICAL  MODE  DISTORTION 

A.  BACKGROUND 

To  continue  the  investigation  of  mode  distortion,  the  single-pass  three 
dimensional  simulation  used  in  Chapter  V  is  modified  to  calculate  the  transverse 
beam  effects  over  multiple  passes.  This  enables  further  research  into  the 
interference  effects  shown  in  Fig  5-8,  and  the  relationship  between  gain  and 
electron  beam  size  as  shown  in  Fig  5-7.  The  multiple  pass  simulations 
introduce  optical  mode  selection  caused  by  the  finite  size  of  the  resonator 
mirrors.  Through  many  passes,  highly  diffracted  modes  may  be  lost  due  to 
mirror  size  and  placement. 

B.  MULTIPLE  PASS  PROPAGATION  AND  INTERFERENCE 

Investigation  of  the  single-pass  optical  interference  pattern  shown  in  Fig 
5-8  over  multiple  passes  will  be  affected  by  the  number  of  passes,  n,  and  the 
mirror  design.  Figure  6-1  shows  the  three  dimensional  simulation  modified  for 
multiple  passes.  For  this  regime,  the  initial  optical  amplitude,  ao  =  0.05,  for  a 
weak  field  system.  At  either  end  of  the  amplitude  plot,  la(x,n)l,  are  curved 
mirrors,  spaced  x^  =  2  apart,  or  twice  the  undulator  length.  The  separation 
distance  of  a  typical  FEL  may  be  larger,  but  this  choice  allows  ease  of 
computation.   The  dimensionless  mirror  radius  of  curvature,  re,  is  determined  by 


the  mirror  spacing,  x^,  and  the  normalized  Rayleigh  length:  Tg  =  t^/2 -i- 2z|/x 


m- 


The  actual  mirror  radius  of  curvature  is  given  by  Lrg.    The  mirror  losses  are 
described  by  O  =  50,  and  the  mirror  edge-loss  is  k  =  0.001. 
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FEL   WAVEFRONTS     (HOLE-COUPLED    OUTPUT) 
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Figure   6-1:  Three   dimensional   sinnulation  for  multiple  passes  and 
destructive  interference. 
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Edge-losses  occur  when  the  mirror  radius,  r^,  is  not  wide  enough  to  capture 
the  higher-order  optical  modes,  losing  them  to  the  non-reflecting  sides  of  the 
resonator.  The  dimensionless  mirror  radius  is  a  function  of  the  Rayleigh  length, 
mirror  spacing  and  edge  loss,  r„  =[ZQ{^ +z^/4zS)  {^/2\oQ{{^ +K)/K))y'^.  The 
amplitude  plot,  \a{x,n)\  (upper-left),  shows  the  evolution  and  growth  of  the 
optical  mode  through  n  =  16  passes,  as  scaled  along  the  base  of  the  plot.  The 
single  contour  line  shows  that  the  amplitude  begins  to  grow  in  the  last  few 
passes.  The  plot  of  resonator  mode  evolution,  la{x,T)l  (center),  shows  the 
entire  optical  mode  at  pass  n  =  16.  The  undulator  ends  are  represented  by  tick 
marks  at  x  =  0  and  x  =  1  along  the  plot  base. 

Careful  analysis  reveals  that  a  radial  "dimple",  similar  to  that  found  in  the 
single-pass  interference  plot  of  Fig  5-8,  is  present  near  the  right  end  of  the 
optical  mode.  In  the  grey  scale  used,  the  destructive  interference  described  in 
Chapter  V  is  present  as  two  unshaded  rays  that  propagate  from  just  outside  the 
undulator  markings  toward  the  right-side  mirror.  The  final  transverse  amplitude 
plot,  la(x,y)l,  shows  the  amplitude  at  the  last  pass,  n  =  16  at  the  end  of  the 
undulator. 

The  electron  phase  velocity  evolution,  /(v.n),  shows  only  small  modulation 
in  weak  fields.  The  phase-space  plot,  (C.v),  shows  slight  bunching  in  this  weak 
field  regime.  Power  evolution,  P{n),  starts  to  grow  at  approximately  the  same 
point  as  the  amplitude  in  la(x,n)l.  Gain,  G,  reaches  a  steady-state  near  n  =S, 
where  G  ~  55%,  and  single-mode  gain. Go  «  yr  predicts  only  28%. 

Figure  6-2  shows  an  FEL  similar  to  that  of  Fig  6-1 ,  but  in  the  strong  field 
regime,  starting  with  Bq  =  5.0,  and  running  over  n  =  160  passes.  The  amplitude 
evolution  plot,  la(x.n)l,   has  a  much  greater  maximum  of  39. 
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FEL   WAVEFRONTS     (HOLE-COUPLED    OUTPUT) 
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Figure     6-2:     Three     dimensional     simulation     without     destnjctive 
interference. 
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The  final  transverse  amplitude  plot,  \a{x,y)\,  and  the  resonator  plot, 
la(x,x)l,  show  an  optical  mode  with  no  radial  "dimple".  The  phase  evolution 
plot,  f{v,n),  and  phase-space  plot,  (C.v),  show  complex  bunching  due  to  strong 
field  saturation.  Power  evolution,  P(n),  starts  near  where  the  amplitude  plot 
begins  to  grow,  and  continues  to  increase  through  the  n  =  160  passes.  If  this 
simulation  continued  for  more  passes,  the  power  evolution,  P{n),  would  achieve 
a  steady  state  value.  The  gain  plot  shows  an  early  peak,  followed  by  a  quick 
reduction  to  near-zero  steady-state  near  n  ~  50.  It  appears  that  the  interference 
"dimple"  is  not  present  at  saturation  in  the  FEL  oscillator. 

C.  GAIN  AND  SMALL  ELECTRON  BEAMS  OVER  MULTIPLE  PASSES 

To  investigate  the  effects  of  electron  beam  size  on  gain  over  multiple 
passes,  the  product  ya|  is  held  constant,  while  the  electron  beam  radius  is 
decreased.  Actual  FEL  gain,  G,  was  numerically  calculated  over  several 
passes,  while  the  single-mode  gain  is  given  by  Go  =  0.135yT  «  yo|.  In  order  to 
determine  whether  the  single-pass  relationship  of  gain  on  electron  beam  size 
shown  in  Fig  5-7  applies  for  multiple  passes  through  the  resonator,  simulations 
were  njn  to  evaluate  gain.  For  these  simulations,  mirror  loss  are  Q  =  1000,  and 
the  edge  loss  is  k:  =  0.001.  The  Rayleigh  length,  zq,  edge  loss,  k,  and  mirror 
spacing,  x^,  result  in  mirror  radius,  r^  =  2.7,  and  radius  of  curvature  Tc  =  0.9. 

Figure  6-3  shows  the  ratio  of  the  actual  gain  to  the  single-mode  gain, 
G/Gq,  plotted  against  the  ratio  of  dimensionless  beam  radius  to  mode  waist, 
Cg/wo  while  holding  jcg  =  0.4.  Fig  6-4  shows  a  similar  plot  for  the  value 
jog  =  0.6.  To  prevent  the  influence  of  high-gain  effects,  the  beam  radius  is 
limited  by  Cg/WQ  >  0.1. 
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Figure  6-3:  Multiple-pass  gain  curve  for  /a|  =  0.4. 
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Figure  6-4:  Multiple-pass  gain  curve  for  ya|  =  0.6. 
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The  curves  resemble  the  single-pass  results  in  Fig  5-7,  with  the  gain 
ratios  decreasing  steadily  as  the  electron  beam  radii  widen  toward  the  size  of 
the  mode  waist.  The  single-mode  gain  product  yo|  is  greater  for  Fig  6-4,  so  the 
gain  curve  for  Fig  6-4  is  higher  than  for  Fig  6-3.  Similar  to  Fig  5-7,  there  is  no 
range  of  values  where  single-mode  theory  is  correct,  only  single  points  where 
G/Go=  1. 
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VII.  CONCLUSIONS 

The  development  of  the  FEL  as  a  useful  tool  continues  through  the 
progress  of  theoretical  and  experimental  efforts  in  several  countries.  The 
potential  for  military  application  of  the  FEL  are  more  promising  than  other 
sources  of  directed  energy.  A  major  challenge  to  all  high  energy  light  sources 
is  the  complex  effects  of  atmospheric  propagation.  Innovative  beam  relay 
methods  may  enable  these  obstacles  to  be  minimized. 

The  proposed  design  of  the  CEBAF  IR  FEL  will  make  important 
contributions  toward  realizing  the  goal  of  a  high  power,  high  efficiency  source  of 
coherent  radiation  for  applications  in  medicine,  industry,  and  possibly  the 
military.  Simulations  of  the  CEBAF  IR  FEL  in  Chapters  IV  and  V  have 
contributed  to  the  FEL  design  selection  process.  The  broad  tuning  range  of  the 
CEBAF  IR  FEL  was  analyzed  for  system  gain  and  power,  and  the  influence  of 
short  electron  pulse  effects.  The  small  electron  beam  of  this  system  created 
strong  distortion  effects  which  were  simulated  in  Chapter  V.  The  interaction  of 
the  optical  mode  and  the  intense  beam  current  resulted  in  higher  FEL  gain  than 
predicted  by  single-mode  gain  theory.  There  exists  a  complex  destructive 
interference  in  the  optical  mode  that  results  in  a  radial  "dimple"  in  the 
propagated  wavefront.  This  phenomenon,  and  system  gain,  was  explored  in 
detail  for  single  pass  simulations  in  Chapter  V  and  multiple  passes  in  Chapter 
VL 
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